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Abstract 

We review the subject of central exclusive particle production at high energy hadron colliders. 
In particular we consider reactions of the type A + B — ^ A + X + B, where X is a fully specified 
system of particles that is well separated in rapidity from the outgoing beam particles. We 
focus on the case where the colliding particles are strongly interacting and mainly they will be 
protons (or antiprotons) as at the ISR, SppS, Tevatron and LHC. The data are surveyed and 
placed within the context of theoretical developments. 
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1 Introduction 



Central exclusive particle production (CEP) is defined as the class of reactions A + B-^A + X + B, 
where the colliding particles A and B emerge intact and a produced state, X, is fully measured. We 
will review the particularly interesting case where the centre-of-mass energy ^/s is large, such that 
one can have rapidity regions Ay > 3 completely devoid of particles between X and the outgoing 
particles A and B (rapidity gaps). With this definition, only high energy colliders are relevant, and 
we shall focus on hadron colliders. The field of central production in nucleus-nucleus collisions at 
the Relativistic Heavy Ion Collider (RHIC) and at the Large Hadron Collider (LHC) merits a review 
in itself and wc will not attempt to cover it. The first hadron collider was the CERN Intersecting 
Storage Rings (ISR), which carried out pp and pp collisions (also dd and aa). The SppS collider 
at CERN followed and, most recently, the Tevatron at Fermilab has examined pp collisions. In 
the near future, the LHC will furnish us with pp (and AA) collisions with much higher energy and 
luminosity than before. It opens up the exciting new possibility of producing heavy systems X, 
such as Higgs bosons, weak vector bosons, supersymmetric (SUSY) particles and other "Beyond the 
Standard Model" (BSM) particles. 

Much of the renewed interest in central exclusive production comes from the higher mass 
reach that will be opened up at the LHC. Added to that are strong constraints that are unique 
to exclusive processes. In some cases these allow mcasiircmcnts of particle properties that can be 
obtained no other way. The "strong constraints" are both kinematic and dynamic. In an exclu- 
sive process such as p + p — )• p + X + p, four-momentum conservation means that knowing the 
incoming proton momenta, and measuring the outgoing proton momenta (in high precision forward 
spectrometers), the mass Mx is determined [1]. This is expected to be possible with a resolution 
of a{Mx) ~ 2 GeV/c^ per event [2], independent of the decay of X, even if X decays with a large 
amount of missing energy, e.g. X W~^W~ — )• e^fi~i'9. With a precise calibration of the forward 
spectrometers (using the QED process — )• p + fi'^fj.~ +p), the central mass can be found with a 

resolution improving as a{Mx)/VN, and if (in the case of a resonance) the width Tx exceeds a few 
GeV, that too can be measured. In addition, the transverse momenta, px, of the scattered protons 
are typically pt ^ 1 GeV/c and so priX) is small. In this case, dynamic constraints, discussed later, 
allow one to also determine the spin J, parity P and charge-parity C of any exclusively produced 
particle, such as a Higgs boson. 

The possibility to instrument the LHC in order to measure the small angle protons produced 
in CEP has been explored through the FP420 R&D project (so named because the forward protons 
would be detected 420 m from the collision region)^ involving a consortium of ATLAS and CMS 
physicists and theorists [2] . Suitable sub-detectors have now been proposed: AFP (ATLAS Forward 
Protons) in ATLAS and HPS (High Precision Spectrometers) in CMS. Precision tracking and timing 
detectors at 220/240 m and 420 m will measure the momenta of both protons. 

In this review our attention will focus on CEP mainly through IP IP — > X and 77 — !■ X, 
although photoproduction, i.e. jlP — >■ X, is briefly mentioned. We use the symbol "iP" to signify 
pomeron exchange (see Section 2). We will neglect contributions from odderon exchange and consider 
data that are not much contaminated by sub-leading Regge exchanges. In classifying the production 
this way we are making explicit the theoretical property that in both QED and in Regge theory the 
scattering amplitudes can often be factorized so as to separate the dependence on the beam particles 

^ There are also plans to install detectors at 220/240 m. 
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from the dependence on the production of the central system. However, it is not always the case 
that the scattering amplitudes factorize this way. In particular, the QCD production of high-mass, 

short-distance (Q^ ^ ^%cd) systems can be computed in perturbation theory (at least up to soft 
rescattering corrections) and the amplitudes do not factorize in quite this manner. Generally, they do 
not factorize wherever the dominance of a single pomeron pole docs not hold. Nevertheless, we shall 
still follow popular nomenclature and refer to the production of a central system via strong dynamics 
as IP IP X, but it is to be understood that this is more of a placeholder than a description of the 
underlying strong dynamics. It is our aim in the following sections to clarify these matters and to 
present the theory underpinning central exclusive particle production. 

In the next section we introduce the relevant ideas in Rcgge theory and show how one 
may describe CEP in terms of pomeron exchange in the f-channel. References [3] and [4] are 
treatises on diffraction and the pomeron and our treatment here will be brief. In Section 3 we 
explain how one can use QCD perturbation theory to compute the CEP of high mass systems. 
While pcrturbative Quantum Chromodynamics (QCD) is the well established theory of high-Q^ 
(short distance) strong interaction processes, diffractive processes almost always involve some low- 

physics, where perturbation theory cannot be used and factorization theorems do not hold. We 
are thus still far from a complete understanding, and that makes for an interesting field, both for 
theory and for experiment. 

We shall then move on to a survey of the experimental data, covering the published results 
from hadron colliders up to the end of 2009. We do not attempt to review CEP in e~^e~ or ep 
machines, and only briefly mention fixed target experiments. We end our review with a forward 
look to the LHC era. 

Some conventions need stating. We will frequently use the term proton to mean p or p 
whenever the distinction is unimportant or obvious. In a reported cross section, such as 2.09±0.90± 
0.19 pb, the first uncertainty is statistical and the second is systematic. Feynman-x, xp = Pz/Pbeam, 
represents the ratio of longitudinal momentum to the beam momentum and ^ = 1 — xp denotes 
the fractional momentum loss of a beam particle (typically a, p or p). The rapidity of a particle is 

defined as y = gin ^ E-p' ) ' pseudorapidity r] = — ln(tan|) with 9 the polar angle; rj = y ioi 

massless particles. 

2 The Pomeron 

A detailed treatment of Regge theory is beyond the scope of this review and we refer instead to 
Refs. [3, 4, 5, 6, 7]. Here we provide a very basic introduction. Pomeron phenomenology was 
developed originally to describe the behaviour of total hadronic cross sections and difFractive/elastic 
scattering at high energies. Most recently, it has been taken to a new level by experiments conducted 
mainly at the HERA ep collider (see Refs. [8, 9] and references therein) and the Tevatron pp collider 
(see Section 7). 

When the CERN ISR came into operation in 1971, the quark model had been formulated 
but not yet QCD [10]; asymptotic freedom was uncovered in 1973. At that time, strong interaction 
physics tended to focus on total cross sections, elastic scattering and few-body reactions (e.g. tt~ + 
p ^ 7r° + n). General methods, based on the unitarity, analyticity and crossing symmetry of the 
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scattering matrix provided the foundations, with Regge theory providing the principal tooL With 
the advent of QCD the emphasis shifted to the investigation of scattering processes at short distances 
for which the strong couphng is small and perturbative methods can be exploited. However, soft 
diffraction and elastic scattering processes cannot be described by perturbative QCD, and Regge 
theory remains an important tool. In Regge theory these process are described as the t-channel 
exchange of "reggeons" {M), which correspond to a sum of mesons (p°,a;", etc.) with the same 
quantum numbers. The contribution of the reggeons to the elastic scattering cross section falls with 
increasing centre-of-mass energy as ^ ^/y/s: where ajpiXt) is the reggeon trajectory which 

is a function of the Mandelstam four-momentum transfer squared, t. By the Optical Theorem, the 
reggeon contribution to total cross sections likewise falls as the centre-of-mass energy increases. The 
observed rise of total hadronic cross sections therefore mandated the emergence of a new reggeon, 
with intercept a]p{0) > 1.0. To generate a non-falling total cross section, the exchange must have 
isospin zero and even charge parity, C = -1-1, i.e. it has the quantum numbers of the vacuum. The 
new reggeon was dubbed the pomeron {IP) after Pomeranchuk, who had previously studied the 
behaviour of vacuum exchange in Regge theory. 

— ? A 

i t 

B • ' B 

Figure 1: Elastic scattering between two hadrons A and B, at centre-of-mass energy ^/s. The 
four-momentum transfer squared is t. 




At high enough centre-of-mass energy, if one assumes the dominance of a single Regge pole, 
the elastic scattering of strongly interacting particles may be described by pomeron exchange, see 
Fig. 1. The elastic scattering amplitude for AB — >■ AB is thus approximated by 



A{s,t) „ „ ... / s 



ajp(t)-l 

S \SoJ 

where 

(2.2) 



ri{t) = i~coi I 1 



is the signature factor, a]p{t) is the pomeron trajectory, iSA.sit) fixes the coupling of the pomeron 
to the external particles and sq is a constant. The Optical Theorem then relates the total cross 
section for AB — > X, (7^, to the imaginary part of the forward {t = 0) scattering amplitude via 

linA{s,0) = s axis) (2.3) 

and so 

/ \ Q2p(0)-1 

aris) ^ f3A{0)f5BiO) i-j . (2.4) 



4 



Thus we see that if the pomeron intercept aplO) > 1 the total cross section rises with energy, in 
accord with the data. Conversely, the contributions of any Regge poles with ajR,(0) < 1 (such as 
those containing the p and tt) become negligible at sufficiently high energy. 

One might hope that the properties of Regge poles, i.e. their intercepts and couplings, 
would emerge from calculations based on QCD. To an extent that is what happens. For example, 
the gluon is known to "reggeize" to a simple Regge pole in perturbative QCD after re-summing to 
all orders in ln(s) (see for example Ref. [3] and references therein). However, calculations are 
generally plagued by the need to focus on processes at short distances, where perturbation theory 
is valid, and total hadronic cross sections certainly do not fall into that class. It is also far from 
clear that amplitudes are dominated by a single Regge pole at high energies, although there is some 
indication that this is so in the case of hadron-hadron elastic scattering at small (but not too small) 
values oit [11, 12]. In that case, fits to data suggest the existence of a pomeron with intercept 

ap{t) « apiO) + a'ptK 1.08 + (0.25 GeV"^) t. (2.5) 

More recent analyses suggest that a global fit to all soft data from the ISR, SppS and Tcvatron may 
require a pomeron with a higher intercept and substantial screening corrections (see for example 
Refs. [13, 14, 15] and references therein). 
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Figure 2: (a) Diffractive excitation of particle i? to a state of mass Mx by pomeron exchange, (b) 
The corresponding cut diagram in the limit of large Mx ■ 



Regge theory is not restricted to the consideration of elastic scattering amplitudes. A 
reggeon calculus can be developed, and used to tackle processes such as those illustrated in Fig. 2(a) 
and Fig. 3(a). Again the dotted lines represent pomerons, and pomeron dominance is presumed to 
pertain if the relevant sub-energies are large enough, i.e. Ay ax ^ 3 in Fig. 2(a) and Ay ax, Aysx ^ 3 
in Fig. 3(a), where AyAx is the rapidity interval between A and X (and similarly for Aysx)- For 
the single diffractive dissociation process represented in Fig. 2(a) one may write [16, 17] 

^x^^^^-^A{t?W)f (^^j OBp{Ml,t) (2.6) 

and we are invited to think of (TBp(M^,t) as the total cross section for BIP scattering at energy 
Mx- It is to be noted that the normalization of aBp{Mjr,t) is a matter of convention. Provided 
Mx is sufficiently large, we expect that it is itself driven by pomeron exchange and aBp{M^,t) oc 
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(M^^)"^'^''^ ^. This is shown in Fig. 2(b) which illustrates the discontinuity in the relevant 
three-body amplitude. We should stress that the pomeron is not a real particle and pomeron- 
induced cross sections are not directly measurable; however they are useful constructs. Going one 
step further, we can rewrite Eq. (2.6) as 



_da_ 
dtd^ 



= fp/A{Lt) aBp(Ml,t) 



(2.7) 



where we define a pomeron "flux" 



flP/A{£.,t)=PA{tfHt)[^ 



2Qjp(t)-l 



(2.8) 



and ^ is the fractional energy lost by the beam particle A, i.e. = ^s. This approach describes 
very well the HERA data on single diffraction dissociation, albeit with a pomeron trajectory that 
differs from that in Eq. (2.5). In particular the t-dependence is consistent with a flat trajectory [9]. 
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Figure 3: (a) Diagram for double pomeron exchange, (b) The corresponding cut diagram in the 
limit of large Mx ■ 



The study of double pomeron exchange (DIPE), illustrated in Fig. 3, has a long history [6, 
7, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27]. In the Regge framework, such exchanges are responsible for 
the CEP process, and we may write the cross section ioi A + B ^ A + X + B in terms of the total 
cross section for two pomerons to fuse, producing the central system X, (ipp: 



da 



diidia dCi dCs 



— flP/AiS,l,ti)fjp/Bi^2,t2) CTPP 



(2.9) 



Again, the are the fractional energy losses, and kinematics fixes M'j^ 
this formula for DPE in Sections 4-6. 



^i^2S- We shall return to 



Eq. (2.7) and Eq. (2.9) clearly exhibit Regge factorization and the similarity to the two- 
photon production case is striking - the pomeron flux playing the role of the Weiszacker- Williams 
flirx in the case of photons (e.g. see Ref. [28]). Unfortunately Regge theory does not tell us how to 
compute the cross section aipip{Mx ,ti,t2), although it does predict the behaviour for large Mx- 
So, although we have a model for the rapidity dependence of the central system, we are not able 
to predict the overall production rate without further model dependence. Furthermore, we should 
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bear in mind that there is no a priori reason why more sophisticated pomeron diagrams should not 
be relevant. 



It is time now to turn our attention to a rather different approach to CEP, namely an 
approach based firmly within perturbative QCD. As we shall see, it does not give rise to Regge 
factorization and leads to essentially different predictions for the CEP of systems for which Mx ^ 
Aqcd- Furthermore, it is an approach that has recently had some striking success in predicting the 
CEP rate of XcO and dijets at the Tevatron. 



3 QCD Models of Central Production 



^1 



Q 




^2 



Figure 4: The relevant lowest order Feynman diagram ioi qq ^ q + H + q. 



3.1 Perturbative QCD 



The perturbative approach to CEP was originally developed in Refs. [29, 30, 31, 32, 33]. We shall 
refer to the approach presented in those papers as the Durham model, and it is our goal in this 
subsection to review the calculation and comment upon its uncertainties. Our present focus is on 
the central exclusive production of a Higgs, i.e. pp — >■ p+H+p, although the theory is not essentially 
different for other high mass systems provided they contain a suitable hard scale. In particular, we 
shall use the same framework when it comes to comparing to Tevatron data on CEP of dijets, di- 
photons and Xc mesons. The calculation starts by considering the process qq ^ q + H + q a,t lowest 
order in perturbation theory, as shown in Fig. 4. The Higgs is produced via a top quark loop and a 
minimum of two gluons needs to be exchanged in order that no colour be transferred between the 
incoming and outgoing quarks. The outgoing quark momenta may be parametrised in terms of the 
momentum fraction each transfers to the Higgs, a;,, and their transverse momenta, q'^j^: 



In 



In 
12 



{l-xi)q'^ + 



(1 - 2:2)92 



(1 - xi)s 

(1 - X2)S 



In 
QlT 



IfJ. 
12T 



(3.10) 



(3.11) 



with s denoting the centre-of-mass energy squared. In the high energy limit the colour- singlet 
amplitude is dominated by its imaginary part, as expected from arguments based on Regge theory 
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(see for example Ref. [3]). It may therefore be computed by considering only the cut diagram of 
Fig. 4. Thus, for small Xi and with the colour-singlet exchange contribution projected out, the 
amplitude is given by 

In the Standard Model, the Higgs production vertex is 

where V = rnj^as/ {4:Trv)F{m'jj/mf), v is the Higgs vacuum expectation value and F ss 2/3 provided 
the Higgs is not too heavy [34, 35, 36]. The Durham group also include a next-to-leading order 
K-factor correction to this vertex [31], though this is extracted from calculations of inclusive Higgs 
production [37, 38] and not from an explicit calculation of the next-to-leading order corrections in 
CEP. 

To proceed further, we parametrise the loop momentum in terms of Sudakov variables, Q = 
aqi + I3q2 + Qt', the (5-functions which fix the cut quark lines on-shell then set a « — /3 « Q|,/s <C 1 
and « = —Q,^. As always, we neglect terms that are suppressed in s, such as the product 
a/3. 



We can compute the contraction V^^^g2 either directly or by utilising gauge invariance. 



which requires that k^Vf^^ = kl^V^I^ = 0. Writing^ hi = Xiqi + kix yields 



1l^tJ.u'l2 ~ ^ ^ >f ~ ^2 '^1T'^2T^hv l-J-J-^J 

since 2ki ■ k2 ~ xiX2-s ~ mfj. Note that it is as if the gluons which fuse to produce the Higgs 
are transversely polarized, i.e. we have for their polarisation vectors: Cj ~ fcjT- Moreover, in the 
limiting case that the outgoing quarks carry no transverse momentum Qt = —k\T = k2T and so 
ei = —£2- This is an important result; it generalizes to the statcmient that the centrally produced 
system must have a vanishing z-component of angular momentum in the limit that the protons 
scatter through zero angle (i.e. q'^ <^ Q\) [29]. Since we are interested in very small angle 
scattering this selection rule is effective. One immediate consequence is that the Higgs decay to b- 
quarks may now be detectable. This is because, for massless quarks, the lowest order qq background 
vanishes identically (it does not vanish at next-to- leading order). The leading order cxchisivc hb 
dijet background is therefore suppressed by a factor ~ ml/mj^. The dominant background thus 
becomes CEP production of ^g-dijets, which can be reduced by 6-tagging both jets. 

Returning to the task in hand, the differential cross section is given by 

klT • k2T 



da fNf.-l\ a° Gp 

X 



(3.15) 



^We can do this because Xi ~ mH/\/s whilst the other Sudakov components are ~ Q^/s 
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with y the rapidity of the Higgs. We are mainly interested in the forward scattering hmit whence 



kiT • k: 



2T 



O^k^ k^ 

IT 2T 



Q4 



(3.16) 



As it stands, the integral over diverges. However, this is because the lowest perturbative order 
is insufficient for this process. Not only must wc include a convolution with non-perturbative parton 
distribution functions (PDFs) but, due to the exclusive nature of the process, the perturbative 
contribution to the amplitude is enhanced at each order in by large logarithms in the ratio 
"^h/Qt- These terms are due to virtual corrections involving soft gluons or partons that travel 
coUinear to the incoming hadrons. They must be summed to all orders in order to give a reliable 
prediction. 



x', Qj 



o 
o 
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o 
o 
o 
o 
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x',Qi 




fg{x,x',Qr,kT,mH) 



Figure 5: The recipe for replacing the quark line (left) by a proton line (right). 



The Durham prescription to account for these effects is shown in Fig. 5. They replace 
ctgCp/T^, at the quark level, with a skewed, imintegrated, gluon PDF, /^(a;, x', Qt, ky, mi/). The 
amplitude is dominated by the region x :s> x' , » k|, and in this limit one may relate the fg 
distribution to the standard, integrated, gluon distribution, g{x,Q,^) [31, 39]: 

fg{x,x',QT,]^T,mH) ^ e-'^^'^/^Rg-^^^ (^^/T{QT,mH)xg{x,ClT)) • (3-17) 

The kr dependence of the scattered proton is assumed to follow a Gaussian distribution, with slope 
parameter, 6 w 4 GeV~^, fixed by a fit to soft hadronic data [40] and cross-checked against the 
ky-dependence of diffractive J/ip production at HERA [41] . The factor Rg is given by 



Rg — 



-^3(2' 2' Qt) 

xg{x, Ql) 



(3.18) 



and accounts for the skewed effect {x ^ x'). Here Hg is the skewed gluon distribution (see for 
example Ref. [42]). In [43] Rg was shown to be given approximately by 



Ro 



22A+3 r(A + 5/2) 
v/^ r(A + 4) 



(3.19) 



if one assumes a simple power-law behaviour for the gluon density, xg{x,Q,'^) x~^. For the 
production of a 120 GeV/c^ Higgs at the LHC, Rg « 1.2 [30]; the off-diagonality therefore provides 
an enhancement of (1.2)'' « 2 to the cross section. Clearly the current lack of knowledge of the 
off-diagonal gluon is one source of uncertainty in the calculation. 
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Equation (3.17) also includes a Sudakov suppression factor, T{Q,T,mH)- This is present 

because real emission from the process is forbidden^. It is the Sudakov factor that collects terms 
in the perturbation series enhanced by logarithms in mj^/Q'^. Taking into accoimt the leading and 
next-to-leading logarithms, i.e. all terms of order a" ln™(m|^/Q|n) with m = 2n, 2n — 1, the Sudakov 
factor has the form [31, 44, 45] 

T{QT,mH)=expi-J^^ -^^^ dz[zPM + nfPM]\ (3.20) 

where are the usual DGLAP splitting functions [46, 47, 48]. It is worth noting that, in order for 
this expression to collect the next-to-leading logarithms, the lower limit on the q^-integral and the 
upper limit on the 2;-integral must be specified precisely. The upper q|, limit in contrast, may be 
multiplied by any 0(1) number without changing the leading or next-to-leading logarithms. This is 
because it corresponds to hard, non-collinear, emission and as such is not enhanced by a logarithm. 
The lower limit on the q|n-integral can be understood as follows. For hard coUinear emissions, 
1 — ^ ~ £)(1), the region < Q|, is already included in the PDFs (see the discussion below). On 
the other hand, for the soft region, 1 — z 0{\qj'\/mH) (relevant for the 1/(1 — z) piece of Pgg), 
emissions with < cannot resolve the exchanged colour-singlet system (the size of which is of 
order I/IQt^I) and so do not contribute. 

The form of Eq. (3.20) is somewhat different from that appearing in much of the literature. 
In particular, in Ref. [33] the upper limit on the z-integral was determined to be |qT|/(|qT|+0.62m//). 
It was shown recently that this determination is incorrect and Eq. (3.20) is in fact the correct form 
of the Sudakov factor [45]. At the present time, the extent to which this modification alters previous 
predictions has not been fully assessed. However, it has been estimated to give approximately 
a factor two reduction in the cross section for Higgs masses in the range 100-500 GeV/c^, with 
the suppression increasing(decreasing) for larger (smaller) Higgs masses [45]. As we shall see, this 
uncertainty still lies within the range of the other uncertainties in the calculation. 

The Sudakov factor solves the Qt — divergence problem we encountered in the lowest 
order calculation; the cross section, integrated over final state hadron momenta, is now given by 

I - 25^^ (/ '-Ws^ (^-(-.«)) (vT«.fe,Q^)))' (3.21) 

and since the exponential Sudakov factor vanishes faster than any power of Qt as Qt ~^ 0, this 
integral is finite. Furthermore, for a 120 GeV/c^ Higgs, the typical Qt ~ 2 GeV/c [49], justifying, a 
posteriori, the use of perturbation theory. Note that the parton momentum fractions are, to a very 
good approximation, equal to the fractional momentum losses of the incident protons, i.e. Xi w 

and yx = ^ln($i/^2)- 

Returning to Eq. (3.17), the derivative structure can be understood as follows (see Ref. [45] 
for more details). Transverse momentum ordered ladder diagrams, as shown schematically in Fig. 6, 
are enhanced^ by logarithms of Q|n. Summing up these contributions evolves the parton distributions 



^Notc that, for a colour-singlct central system, a single real emission is forbidden regardless of experimental cuts, 
since it violates colour conservation. 

*In a physical gauge the large logarithmic corrections are organised into ladders on a diagram-by-diagram basis. 
In contrast, the ladder structure only emerges in a covariant gauge after a sum over diagrams and application of the 
Ward identity (see for example Ref. [50]). 
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associated with the upper and lower hadrons to the scale Q^. The central hard process however, is 
sensitive to the transverse momenta of the final rung in each ladder. If the final rung has an energy 
much larger than its transverse momentum, E ^ |Qt|, then it is included in the PDF and so we 
require the gluon PDF unintegrated over transverse momentum, dg{x,Q,^)/dCl^. If, on the other 
hand, the final rung is soft, E ^ |Qt|, then it may not be included in the PDF. This contribution 
is accounted for by the derivative of the Sudakov factor. Shown also in Fig. 6 (by the blob labelled 
T) are the corrections making up the Sudakov factor. These are virtual corrections to the gg — ^ H 
sub-process with larger transverse momenta than the final ladder rungs. 






Figure 6: Schematic diagram of the transverse momentum ordered ladder corrections which evolve 
the PDFs to the scale Q^. Solid lines denote either quarks or gluons. 



So far we have only discussed the Higgs production case, however the general features of the 
calculation are expected to remain the same for other central systems. The differential cross section 

for the production of a central system X , of invariant mass v^, may therefore be written as 



da 



dC 



(3.22) 



The partonic cross section, a, is related to the matrix element for two on-shell gluons to produce 
the central system as 



da{gg X) ^ j,\M{gg ^ X)\' dPSx 



(3.23) 



where dPSx is the phase-space of the final-state, X and 

1 1 

aia2 A1A2 



M{gg^X) 



2N^ 



Y E E ^a,aJx.x,Ml[ll{gg ^ X) 



(3.24) 



with A^^^'J!, the amplitude for two on-shell gluons, with colours and helicities Ai, to fuse to 
produce X. Finally, the effective luminosity is given by 



d£ 

dsdydqfrpdq^2T ^ 



s \ N^- 1 J 'Qr'/g(^i''^i'Qr'qiT:^)/g(^2,4.QT,q2T.^)j • (3-25) 



Returning to Higgs production, before we can compute the cross section we need to intro- 
duce the idea of rapidity gap survival. The Sudakov factor has allowed us to ensure that the exclusive 
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nature of the final state is not spoilt by perturbative emission off the hard process. What about 
non-perturbative particle production? The protons can in principle interact independently of the 
perturbative process discussed above, and this interaction will lead to the production of additional 
particles. We need to account for the probability that such emission does not occur. Provided the 
hard process leading to the production of the Higgs occurs on a short enough timescale, we may sup- 
pose that the physics which generates extra particle production (including additional parton-parton 
interactions) factorizes, and that its effect can be accounted for by an overall factor multiplying the 
cross section we just calculated. This is the "gap survival factor", 5^, introduced in Ref. [51] (see 
also [52]). It is defined by 

da{p + H + p\no soft emission) = d(j{p + H+p)x S"^ , (3.26) 

where da{p + H + p) is the differential cross section computed above. The task is to estimate 

S'^. Clearly this is not straightforward since we cannot utilize QCD perturbation theory. However 
data on a variety of processes observed at HERA, the Tcvatron and the LHC can help improve 
our understanding of rapidity gap survival. We expect that 5*^ ^ 3% for CEP at the LHC (sec for 
example [53] for an overview). Early measurements at the LHC of rapidity gap processes will provide 
valuable information on S"^. Note that the gap survival probability has some process dependence: it 
will be higher for large impact parameter, r, interactions. Thus we expect that S'^(77) > 5^(7JP) > 

The issue of gap survival is frustrating to theorists. The spacetime structure of the collision 
suggests that it may approximately factorize. Possible non-factorizing contributions have been 
examined in Ref. [54], although their importance has been questioned in Refs. [55, 56, 57]. Assuming 
it to factorize, we may use data to constrain it and accurate theoretical predictions are possible. 
Here we shall present a simple model of gap survival which should provide a good starting point for 
understanding the more sophisticated treatments in Refs. [14, 55, 57, 58, 59, 60]. 

Dynamically, one expects that the likelihood of extra particle production will be greater if 
the incoming protons collide at small r. The simplest model which is capable of capturing this feature 
is one which additionally assumes that there is a single soft particle production mechanism, let us call 
it a "rescattering event", and that rescattering events are independent of each other for a collision 
between two protons at transverse separation r. In such a model we can use Poisson statistics to 
model the distribution in the number of rescattering events per proton-proton interaction: 

Pn(r) = ^exp(-x(r)) . (3.27) 

This is the probability of having n rescattering events where x(^) is the mean number of such events 
for proton-proton collisions at transverse separation r. Clearly the important dynamics resides in 
x(r); we expect it to fall monotonically as r increases and to be much smaller than unity for r much 
greater than the QCD radius of the proton. Let us for the moment assume we know x(^)) then we 
can determine S"^ via 

2 _ / dr da{r) exp(-x(r)) 



/ dr da{r) 



(3.28) 



where da{r) is the cross section for the hard process that producc;s the Higgs expressed in terms of 
the transverse separation of the protons. Everything except the r dependence of da cancels when 
computing 5^ and so we need focus only on the dependence of the hard process on the transverse 
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momenta of the scattered protons (q^x), these being Fourier conjugate to the transverse position of 
the protons, i.e. 



da(r) o, [ij d^q'i^ e'i---/2 exp(-6qf^/2)) x ( J d^q^^ e'^'i---/^ exp(-6q'2?r/2))]^ 



exp(^--J . (3.29) 

Notice that b here is the same as that which enters into the denominator of the expression for the 
total rate (Eq. (3.21)). There is thus a reduced sensitivity to b when one includes the soft survival 
factor, due to the fact that as b decreases so does S"^ (since the collisions are necessarily more central). 
Thus what matters is the ratio S'^/b'^. 

It remains for us to determine the mean multiplicity x(^)- If there really is only one type 

of rescattering event^ independent of the hard scattering, then the inelastic scattering cross section 
can be written 

o-ineiastic = j d^r(l - exp(-x(r))), (3.30) 
from which it follows that the elastic and total cross sections are 

cTeiastic = j d^r{l - e^'Xir) (3.31) 
cTtotai = 2 j d2r(l - exp(-x(r)/2)). (3.32) 

There is an abundance of data which we can use to test this model and we can proceed to perform 
a parametric fit to x(r). This is essentially what is done in the literature, sometimes going beyond 
a single channel approach. Suffice to say that this approach works rather well. Moreover, it also 
underpins the models of the underlying event currently implemented in the PYTHIA [61, 62] and 
HERWIG [63, 64] Monte Carlo event generators which have so far been quite successful in describing 
many of the features of the underlying event, see Refs. [65, 66, 67, 68, 69, 70, 71]. Typically, models 

of gap survival predict to be a few percent at the LHC. Although data support the existing 
models of gap survival there is considerable room for improvement in testing them further, thereby 
gaining greater control of what is perhaps the major theoretical uncertainty in the computation of 
exclusive Higgs production. Early diffractive studies at the LHC should give valuable information 
on S'^ (see for example [72] ) and comparisons between different CEP channels at both the Tevatron 
and LHC will further constrain the theory (see Ref. [73]). 



3.2 Other Models 

Bialas and Landshoff (BL) developed [74] a model inspired by a blend of QCD perturbation theory 
and Regge phenomenology to predict the rate of central Higgs production m pp ^ p + H + X + p, 
i.e. central inclusive Higgs production. Although their calculation looks ostensibly as if it is valid 
for exclusive production, BL were careful to emphasise that "additional... interactions... will generate 
extra particles... Thus our calculation really is an inclusive one", but with both protons having 

® Clearly this is not actually the case, but such a "single channel eikonal" model has the benefit of being simple. 
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xp > 0.95. Inspired by this calculation, the Saclay group [75] attempted to compute the rate for 
pp ^ p + H + p under the assumption that the exclusive rate can be obtained, modulo an overall 
gap survival factor, from the calculation in Ref. [74]. We note in passing that the first attempt 
to compute the central inclusive Higgs cross section, for proton-proton collisions, appeared in [76]. 
This was followed by a similar analysis for nucleus-nucleus scattering in [77], where the first effort 
to estimate the exclusive Higgs cross section was made. Other early studies of central inclusive and 
exclusive Higgs production can be found in Refs. [78] and [79] respectively. 

We can understand the BL calculation starting from Eq. (3.15). BL account for the coupling 
to the proton in a very simple manner: They multiply the quark level amplitude by a factor of 9 
(which corresponds to assuming that there are three quarks in each proton able to scatter off each 
other). Exactly like the Durham group they employ a form factor, exp{—bq'^^), for each proton (at 
the cross section level), with 6 = 4 GeV~^. Since BL are not interested in suppressing radiation, 
they do have a potential infrared problem since there is no Sudakov factor. They dealt with this by 
following the earlier efforts of Landshoff and Nachtmann (LN) in replacing the perturbative gluon 
propagators with non-perturbative ones [80, 81]: 



fc2 



^ ■ Aexp(-A;VM^)- (3.33) 



Rather arbitrarily, = An was assumed, except for the coupling of the gluons to the top quark 
loop, where Ug =0.1 was used. 

Following LN, n and A are determined by assuming that the pp elastic scattering cross 
section at high energy can be approximated by the exchange of two of these non-perturbative gluons 
between the 3x3 constituent quarks. The imaginary part of this amplitude determines the total 
cross section for which there are data to which they can fit. To carry out this procedure successfully, 
one needs to recognize that a two-gluon exchange model is never going to yield the gentle rise with 
increasing centre-of-mass energy characteristic of the total cross section. BL therefore also include an 
additional "reggeization" factor of in the elastic scattering amplitude, where the pomeron 

trajectory is given by Eq. (2.5). In this way the two-gluon system models pomeron exchange. They 
found that /i « 1 GcV and A « 30 GcV~^ gave a good fit to the data. Similarly, the amplitude for 
central Higgs production picks up two reggeization factors. 

The inclusive production of a Higgs in association with two final state protons is clearly much 
more infrared sensitive than the exclusive case where the Sudakov factor saves the day. Nevertheless, 
the Saclay model docs not include the Sudakov suppression factor. Instead it relies on the behaviour 
of the non-perturbative gluon propagators to render the Qt integral finite. As a result, the typical 
Qt is much smaller than in the Durham case. 

Putting everything together, the Saclay model of the cross section for pp^ p + H + p gives 



^ ?1 S2 



) J -^:^Qt ( ^ 1 exp(-3QT/M )- 



2 



(3.34) 



The only difference between this and the original BL result is the factor of S^. Integrating over the 
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final state transverse momenta and simplifying a little gives 



9cr 2 

&y ^ h + 2a' ln(l/a) b + 2a' ln(l/^2) 

H A^cW ^/2(27^)M47^)4 Uij 5^ 

A very similar formula was also presented by Bzdak [82], with the important difference that the 
Sudakov factor was reinstated, i.e. the integral over in Eq. (3.34) is replaced as follows 

j d^Qr Qt exp(-3Q|./At2) -)> j d^Qr Q| exp(-3Q|/At^) T(QT,m//) (3.36) 

with T{QT,mH) given in Eq. (3.20). 

Not surprisingly, the Saclay and Bzdak approaches produce very different predictions: the 
Sudakov factor strongly suppresses the value of the predicted cross section. Broadly speaking, the 
Saclay prediction is similar to the Durham one for Higgs masses aroimd 100 GeV/c^ and (due to 
the absence of a Sudakov factor) overshoots it at higher masses. In contrast, the Bzdak calculation 
is typically much smaller than the Durham result. In both models, the choice of an exponentially 
falling gluon propagator means that there is little place for a perturbative component. However, 
as the Durham calculation shows, there does not seem to be any good reason for neglecting con- 
tributions from perturbatively large values of Qt- In addition, neither calculation can lay claim 
to a systematic summation of the leading double and single logarithms, which is in contrast to the 
Durham calculation. As we shall see in Section 7.2, the CEP dijet data from CDF [83] appear to 
exclude the Saclay model in favour of the Durham one. 

The perturbative Sudakov factor is also included in the approach of Refs. [84, 85, 86]. 
This latter approach also uses perturbative gluons throughout the calculation but Regge factors are 
included to determine the coupling of the gluons into the protons, rather than the unintegrated 
partons of the Durham model. The results are broadly consistent with those of the Durham model. 

Finally, we mention the work of Szczurek and Lebiedowicz, which takes the Durham model 
as the basis for (part of) their predictions for pp p + /o(1500) + p [87] and pp —> p + Xc + P [88] . 

In the following sections we shall survey existing results on CEP from experiments performed 
up to the end of 2009. The reader most interested in the LHC may skip Sections 4-6, but Section 7 
remains relevant. 



(3.35) 



4 CERN Intersecting Storage Rings (ISR) 

The ISR started in 1971, and provided pp and pp (also dd and aa) collisions from = 23 GeV to 63 
GcV, greatly extending the energy range above the ^/s = 7.6 GeV of the CERN PS and Brookhaven 
AGS. The higher energies allowed for a more direct study of pomeron exchange. Ganguli and Roy 
gave an excellent review of Regge phenomenology as it stood after the first eight years of ISR 
operation [7] . Before discussing the data, a note on nomenclature is in order. We will, especially in 
this section, follow the historical terminology and use the term "double pomeron exchange" (DIPE) 
for inclusive studies of central particle production or for exclusive particle production when we wish 
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to emphasise its IP + IP X component. At the low masses typical of the ISR, the central system 
will be dominated by production of tt+tt^ making it in effect exclusive. We shall aim to reserve 
"CEP" for strictly exclusive central systems. It is now accepted that the pomeron is predominantly 
gluonic, which means that DIPE should be a good place to look for glueballs, and at the ISR (and 
even at the lower ^/s of the SPS) this became a strong motivation to study such processes. This 
section provides some historical background, covering the ISR experiments that published results on 
DIPE and CEP. 

4.1 Seairches for double pomeron exchange 

At the start of the ISR era, light hadron spectroscopy was being studied intensively and it was 
"known" that the produced particles have low transverse momentum, pT- Feynman had just invented 
the parton model and proposed that Lorentz-invariant inclusive particle production cross sections 
should "scale", i.e. Ed^a /A^p{pttXfjs) should become independent of s. Several single particle 
spectrometers tested this and scaling was found to be approximately true at small polar angles 9, 
but was dramatically broken at > 2 GeV/c due to hard scattering of partons, with the growing 
realisation that these are quarks and gluons, interacting according to the rules of QCD. 

The ISR physics with which we are concerned in this review has two primary aspects: 
understanding the pomeron and looking for glueballs. A full understanding of strong interaction 
physics must necessarily include an understanding of the pomeron/glueball sector. Yet, even after 
almost 40 years of QCD, the physics of these objects is still not well understood. The ISR enabled 
for the first time the study of pomerons beyond simple elastic scattering and low mass diffractive 
excitation. One could study reactions such as JP + p ^ X for Aix > 10 GeV/c^ and, our main 
concern here, IP + IP ^ X. At the ISR, the state X could be separated in rapidity from both 
outgoing coherently scattered protons by large gaps Ay ^ 3. 

In colliding proton beams the total rapidity coverage is given by ln(s/mp), which is 8.4 at the 
top ISR energy of ^/s = 63 GeV. It was found that single diffractive excitation, e.g. p+p — ?• p+p* — )• 
p+ (p7r+7r~) as studied at lower energies, now extended well above the resonance (N*) region. The 
inclusive forward proton spectrum has a distinct peak [89] ioi xp = 1 — ^ > 0.95, independent 
of -^s (Feynman scaling) and corresponding to IP exchange. The diffractive mass calculated from 
Xi? = 1 — M"^ I s extended up to 14 GeV / c^ , with this limit arising from the requirement of a rapidity 
gap adjacent to a leading proton of Ay > 3 units. For CEP, we similarly find that requiring both 
protons to have xp > 0.95 (or equivalently two gaps of Ay > 3) one is restricted to Mx ^ 3 GeV/c^. 

Earlier searches for DIPE in fixed target experiments, such as that by the France-Soviet 
Union (FSU) Collaboration at Serpukov [90] , were only able to put upper limits on the cross section. 
FSU used a bubble chamber to study the exclusive reaction p + p ^ p + tt+tt" + p at Pboam = 69 
GeV/c, i.e. ^/s = 11.5 GeV. All the events were found to have both pions close in rapidity to one 
of the outgoing protons, consistent with single diffractive dissociation. A limit of (Tdpe ^ 20 /ib 
was determined for events with M^+„- ^ 0.7 GeV/c^ and \y^+ +y.,r- 1 < l-6i a factor ~ 20 less than 
the dominant single diffractive cross section (j{p + p ^ p + (ptt+tt")). At this low value of \/i the 
beam rapidity is only ybeam = 2.5, too low for the distinctive kinematics of DIPE. Another bubble 
chamber experiment, at Fermilab [91], with the higher beam energy of pbeam — 205 GeV/c (Vs = 
19.7 GeV), placed an upper limit of 44 fjh on the DIPE cross section. It required the higher ISR 
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energy for the DIPE signal to emerge. 

Chew and Chew [23] made a prediction for the DIPE signal, based on the single diffraction 

dissociation cross sections ioi A + B ^ A-\- X [a^) and A + B ^ B + X (a^) and an assumption 
of factorization, following Mueller [92]. They anticipated a differential A + B^A + X + B cross 
section given by 

da _ 1 da"^ da^ 

dtAdted^Ad^B ~ ariAB) dt^d^A disd^s ' ^' 

where ^a,b are the fractional momentum losses of the colliding particles and ariAB) is the total AB 
cross section. Integrating over tA,B and over the ^a,b range expected to be dominated by DIPE, 
they predicted aD]PE{'^~^'^~) ~ 65 fxh at y/s = 63 GeV, mostly with Mx < 1 GeV/c^. 
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Figure 7: Experimental DIPE cross sections (/ib) versus s (GeV^) at the ISR together with the 
Regge calculations of Ref [26]. The full circles and the rising solid line arc for two gaps with Ay > 3. 
The dashed line is for Ij/ttI < 1-0 and the dotted line for < 1.5. Figure from Ref. [26]. 



A more sophisticated calculation using Eq. (2.9), and illustrated in Fig. 3, was performed 
by Desai et al. [26] (see also Ref. [7]). To estimate app they assumed pomeron dominance, i.e. they 
took 



where 



^,2/^^„ tnj2 , ^ ^„,2 u ^ _ 9P]Plp{tl)gFPlp{t2) f ^aJ,(0)-l 

ap 



(4.38) 
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g]Pipp{t) is the triple-pomeron coupling (e.g. see Fig. 3(b)) and ap is the pomeron contribution to 
the total cross section, i.e. 

/ \ Q!ip(0)-1 / \asi{Q)-l 

aT{AB)=api-\ ^""^WJ ' ^^'^^^ 

The values of gppF{t) and the Regge trajectory parameters were taken from the analysis of Field 
and Fox [93] and secondary [IR) exchanges were included as a correction (around 30-50% at the 
lowest vahies of s, falling to under 10% at the higher values). They calculated cross sections both 
for two gaps with Ay > 3 and for fixed ranges, < 1.0 or \yT^\ < 1.5, of the central pions, through 
the ISR energy range. The predictions are shown together with the ISR data in Fig. 7 and the 
agreement between the two is rather good. The predicted cross section for two gaps of fixed length 
rises with ^/s, while for a fixed l^jrl < Umax cut the cross section falls with ^/s. The predictions also 
indicate that the two protons should each have an exponential t-slope with hupE ^ ^eiastic/2 and 
they should be uncorrelated in azimuth and t. In the ISR energy range ^elastic ~ 13 GeV~^, with 
a small f-dependence, rising like h{s) = c + dln(s). Note that, under the assumption that <jpp is 
independent of Af|-, the Desai et al. approach is almost identical to that of Chew and Chew, the 
only substantive difference being the treatment of secondary exchanges. 

Desai et al. also consider a pion exchange model for upp, and in this context they point out 
that absorptive corrections may well be important. We refer to Ref. [26] and references therein for 
more detail. We shall now turn to discuss the data, and note that Fig. 7 contains the main results 
on the integrated cross section measurements. Table 1 summarises the ISR DIPE cross section 
measurements. 

The first experimental evidence for DIPE, at the ISR, was presented by the ARCGM Col- 
laboration [94] which had forward proton tracking (but no momentum measurement) and full an- 
gular coverage with scintillation counters. Candidate events had two leading non-coUinear proton 
tracks with polar angle 6 < < 10 mrad, and exactly two hits in a scintillation counter ho- 
doscope covering |r;| < 1.5. Approximately 100 events were found at each ISR \/s value. Setting 
tA,B — — (Pbeam^^,B)^ they found no evidence for any correlation between Ia and t^, and writing 
d?a /dtA^tB = a e~''(*'^+*®^ they measured a slope b = (9.9 ± 1.8) GeV~^, not significantly changing 
with energy. This is somewhat higher than 6eiastic/2 but is consistent at the 2a level. The cross sec- 
tion is also approximately constant, from (Tdpe = (28.4± 8.1)/ib at s/s = 31 GeV to (20.2 ± 3.3)/ib 
at ^/s = 63 GeV. These are expected characteristics of DIPE. 

The Split Field Magnet (SFM) facility at the ISR had nearly full angular coverage for 
charged particle tracks using proportional chambers, with dipole fields in the forward directions to 
analyse the momenta of the scattered protons. The central field was complicated (a quadrupole) 
and central particles were not identified but were assumed to be pions. A study by the CCHK 
Collaboration [95] selected events with two leading positive tracks (a;^ > 0.9) and forward rapidity 
gaps Ay > 2. They observed the DIPE characteristics of no correlations in </> and t, and measured 
b = (5.5 ± 0.9) GeV~^. Extrapolating in \t\ they found cjdipe = (25 ± 10)/ub, consistent with 
the subsequent Desai et al. calculation. The mass distribution (assuming the central particles to 
be pions) had mostly M^+^- < 1 GeV/c^ with no significant p tt+tt^ signal, indicating DIPE 
dominance, as IP IP — > p is forbidden by isospin. The authors point to the absence of an /o(980), 
but we will see later that it appears strongly in CEP as an edge, but not a peak. Au et al. [96] 
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later proposed that this is because the tt+tt ^ tt+tt cross section is already as large as allowed 
by unitarity up to M^^r ~ 1 GeV/c^. 

The CHOV Collaboration also studied p + p ^ p + tt+tt" + p at the SFM [97], selecting 
events with xp > 0.9 for both protons and \y^\ < 1.0. After applying a four-constraint kinematic 
fit, the A(j){pp) distribution showed two components, attributed to single diffraction and DIPE. The 
inclusive proton spectra showed a diffractive peak iov xp ^ 0.95, while the region 0.90 < xp < 0.95 
was dominated by non-diffraction (Reggeon M exchange). They claimed, for \y^\ < 1.0, ctdpe = 
(7.1 ± 1.0)/xb at 23 GcV and (6.0 ± 1.5)/(xb at 45 GeV. Assuming the central particles to be pions, 



The CCHK Collaboration then followed the CHOV analysis with a more detailed study 
using five ISR energies and much higher statistics [98]. A sample was obtained with xp > 0.9, 
predominantly M^+^- < 1 GeV/c^, with no observed p resonance and with an indication of a 
dip around 1 GeV/c^. The helicity angle, defined as the angle in the tt+tt" centre-of-mass frame 
between the 7r+ and the direction of the tt+tt^ system in the overall centre-of-mass frame, showed a 
flat distribution which suggests a dominantly J = system. No correlation was seen \n(j) ov t between 
the protons, and the slope was measured to be 6 = 7.0±0.5(6.1±0.8) GeV-=^ at = 30.7(52.8) GeV. 
Extrapolating in \t\, requiring both protons to have Xp > 0.9 and to have both rapidity gaps Ay > 3, 
the cross section (tdipe{'^~^t^~) was found to rise from (3.6±1.7)/ib to (11.7±3.0)/ib as y/s increases 
from 30.7 to 62.3 GeV. 

The CHM Collaboration [99] studied events with two leading protons, measured with single- 
arm magnetic spectrometers with a{p)/p ~ 0.6%, at ~ 30 mrad (and excluding elastic scattering 
by acollinearity) . A scintillation counter hodoscope covered 98% of the full solid angle and measured 
the directions of central charged particles, but without magnetic analysis. The mass and rapidity of 
the central hadronic system were determined entirely from the measured protons via 



However, the study of centrally produced resonances was limited by the lack of tracking and the 
modest mass resolution of (Tmx ~ 200 MeV at ^/s = 30 GeV. 

Table 1 gives an at-a-glance summary of the above experiments; more details are given in 
the text and, of course, in the original publications. 

4.2 Searches for glueballs 

Glueballs (G) are gluon bound states that are widely anticipated in QCD. There is general agreement 
that the lightest glueball should be a scalar with mass in the range 1 — 2 GeV/c^, with pseudoscalar 
and tensor glueballs at higher mass. However, understanding the light scalars has proven to be 
experimentally and theoretically very challenging, and after thirty years of investigation the glueball 
sector is still not well understood, due in no small part to the likely mixing of any scalar glueball 
with qq (meson) bound states of the same quantum numbers, ft is not our aim here to review the 
experimental status of glueballs; excellent recent reviews can be found in Refs. [102, 103]. Rather we 
will focus our attention on data collected in DIPE that is of relevance to the search for glueballs. It 
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Table 1: ISR experiments reporting DIPE cross sections. See text and cited papers for details. The 
forward protons arc always tracked, but their momenta were not always measured. Only the AFS 
experiment identified the central particles, otherwise pions were assumed. "For M^^ < 1 GeV/c^. 
''For 1 GeV/c2 < M^^ < 2.3 GeV/c^. '^Assuming slope 6 = 6 GeV'^. "^For Ay > 3. 
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was Robson [104] who, in 1977, first suggested that DIPE should provide a gluon-dominated channel 
that favours central glueball production. 

Glueball searches in CEP are not only favoured by the ghion-dominated nature of IP, but 
also the quantum numbers the central state must have: zero isospin, CP = ++ and even spin; it is 
thus a quantum number filter. The lightest known states satisfying these constraints are the very 
broad (and long controversial) /o(600) (or a) and the /o(980), which lies at the KK threshold. 
Note that scalar glueballs cannot lie on the pomeron trajectory aip{t) specified by Eq. (2.5), which 
instead anticipates a J = 2 tensor glueball with Mq ~ 2.1 GcV/c^. 

Following Robson's suggestion, Waldi, Schubert and Winter examined thcp+7r+7r^+p data, 
taken by the CHOV Collaboration, for resonant structure [105]. They observed a bump between 
1100 and 1500 MeV/c^, which they incorrectly attributed to the /2(1270). Later the Axial Field 
Spectrometer (AFS) experiment (which we discuss in more detail in the next paragraph) showed 
that this bump, visible in Fig. 8, is predominantly S-wave and is broader than the /2(1270), with 
a small ( ^ 10%) D-wave contribution from the /2(1270) [101]. Since the /2(1270) is also seen 
in e+e~ e+ + tt+tt" + e~ (two-photon) production with a cross section that favours a large qq 
composition, it is not now considered to be a glue-dominated state. 



Exclusive 7t*7t' 




M(7ttt) GeV 

Figure 8: M^+„- spectrum in DIPE at the ISR (Axial Field Spectrometer, R807 [100, 101]). Figure 
from Ref. [101]. 

Compared to the CHOV analysis, the AFS Collaboration benefitted from a factor 20 more 
statistics (89,000 exclusive events), better mass resolution, cjm = 10(25) MeV/c^ at 1(2) GeV/c^, 
and particle identification [100, 101]. Forward proton track detectors and rapidity gap detectors were 

added to the central spectrometer, and CEP processes with X = n'^ 7r~ , K~^K~ ,pp, and 7r"'"7r~7r+7r~ 
were measured. The two-pion mass spectrum M^+^- shown in Fig. 8 shows several interesting 

features. The cross section rises steeply from threshold and there is no sign of p° — > tt+tt", which 
is forbidden in DIPE. There is however an order-of-magnitude sharp drop just below 1 GeV/c^, 
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followed by another bump-dip structure. The absence of a p, which is prominent in lower energy 
p + p — > p + TT+TT" +p experiments (e.g. at the SPS [106]), is a sign that M exchange is unimportant 
in the AFS data. Since photons can be exchanged across large rapidity gaps, they could give rise to 
a centrally photoproduced p via 7-I-JP p. However the cross section is much lower than for DIPE, 
and protons with \t\ < 0.01 GeV^ were not detected, which further diminishes the rate. Despite the 
absence of a p signal in the mass distribution, when the angular distributions were studied and a 
P-wave (J = 1) component was selected, a small p signal was seen, demonstrating the sensitivity 
of such distributions^. Conversely when the fixed-target experiment WA102 [106] projected out the 
S-wave their p peak disappeared and the M„+„- spectrum at y/s = 29 GeV had a similar shape 
to the AFS spectrum. The partial wave analysis of the AFS group showed that the ( J=0) S-wave 
dominates as far up as 1700 MeV/c^ [101]- 

The striking drop at M^+^- « 1 GeV/c^ in Fig. 8 is the /o(980), the lowest mass narrow 
scalar state with full width 40 — 100 MeV/c^ [107]. According to Au et al. [96] the interaction can 
be understood as iP + iP — )■ tttt — tt+tt" with the final state interaction phase shift dominating the 
structure. Coupled channels such as intermediate KK states contribute; similarly in the K^K~ 
final state intermediate tttt states IP + JP ^ tttt/KK — )■ K~^K~ need to be taken into account. Au et 
al. claimed that the coiiplcd channel analysis could be interpreted as showing two narrow states near 
1 GeY/c^: a KK "molecule" and a glueball candidate. Narison [108] suggested from an analysis 
of widths and couplings that both /o(600) and /o(980) have about 50% (gg) and 50% -^{uu + dd) 
in their wave functions. Mennesier, Minkowski, Narison and Ochs [109] analysed tt+tt" scattering 
below 700 MeV/c^ in an analytic K-matrix model. This is the region of the broad /o(600). They 
argue that the small direct coupling of this state to 77 is "hidden" by the tttt, i.e. /o(600) — >■ tttt — >■ 77 
and find that while r(/o(600) — > 77) ^ 3.9 ± 0.6 keV, large enough to suggest a qq nature, after 
taking into account the rescattering, the direct width is only r(/o(600) — )■ 77) ~ 0.13 ± 0.05 keV. 
This is considerably smaller than the width anticipated in the non-relativistic quark model for a 
{uu + dd) meson of similar mass. However, the interpretation of this state as a glueball is far 
from clear. For example, it is also commonly identified as the scalar partner of the pions under 
spontaneously broken chiral symmetry, especially in the context of NJL models [110]. 

Minkowski and Ochs in particular have interpreted the scalar tttt data from this and other 
experiments in terms of a single very broad state extending from 400 MeV/c^ to about 1700 MeV/c^, 
with the /o(980) and /o(1500) states superimposed, interfering destructively and therefore manifest- 
ing as dips in Fig. 8. They find that in this energy range, after the /o(980) and /o(1500) are 
subtracted, the tttt clastic amplitude describes a full loop in the Argand diagram. They take this 
broad object, which they call the "red dragon", to be the lightest scalar glueball [111]. This is a 
remarkable claim: The full tt+tt^ spectrum of Fig. 8 from threshold to ~ 1700 MeV/c^ is one very 
broad scalar glueball. However any state with a width F exceeding a few hundred MeV decays 
within ^ hc/T < 1 fm, and will not propagate beyond the interaction region as a free particle. The 
situation remains unclear however: Spanier, Tornqvist and Amslcr reviewed the scalar mesons in 
the PDG [107] and concluded that "The 1 = 0, J^*-^ = 0++ sector is the most complex one, both 
theoretically and experimentally." There is however something of a consensus that the scalar meson 

^This method could in principle be applied to, e.g., exclusive p + p — > p + W'^W~ + p (or X = 66) at the LHC, 
if one had high enough statistics and good angular coverage. A scalar state (e.g. a Higgs) would show up only in the 
S-wave component. 
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sector does contain a scalar glueball degree of freedom; the challenge is to understand its mixing 
with the quark states [112]. 

The AFS experiment also measured central K'^K~ production [100, 101]. The cross section 
rises very steeply from threshold, probably due to /o(980), but with no sign of exclusive (p — > , 
which is forbidden in DIPE. The cross section ratio a{K'^K~)/(T(Tr^Tr^) even exceeds 1.0 above 
threshold, perhaps an indication of a gluonic component. Exclusive p + (pp) +p with mass Mpp from 
threshold to 2.8 GeV/c^ was also seen (64 events), but with no evident structures. An extrapolation 
in \t\ gives an estimate of the total a{p + p p + pp + p) « 35 nb (see Table 1). Wc are not aware 
of a prediction for exclusive baryon pairs. One could expect a similar cross section, but scaled by 
M~^, for other baryon pairs, e.g. AA, as the pomeron is flavour-blind. DIPE production of hyperon 
pairs was not studied at the ISR, but it would be possible at the Tevatron and the LHC. 

The ISR also provided aa collisions with ^/s = 126 GeV. The forward drift chambers 
installed for the CEP measurements could distinguish charge Q = 1 and Q = 2 particles. The 

exclusive process a + a — > a + 7r"'"7r~ + a was seen [101], with the spectrum having the 

same shape as in pp collisions within limited statistics (395 events) . The a-particles must have been 
coherently scattered, making this is a very clean channel for DIPE, but the statistics were too low 
for serious spectroscopy. 

The ISR was terminated in 1983, after the SppS collider was turned on. From a CEP per- 
spective, what had we learnt? Well, the existence of double pomeron exchange had been established 
and was in agreement with Regge predictions, but it had been a long struggle. While CEP with 
large rapidity gaps in hadron-hadron collisions can come from 77,75-*, or IPJP interactions, only 
the latter had been observed. There was no evidence for the odderon, the C = — 1 partner to the 

IP, except for a small difference in pp and pp elastic scattering [113] at t = —1.2 GeV^ (the dip 
region). The reach of CEP was up to Mx » 3 GeV/c^, which was favorable for hadron spectroscopy 
and especially for glueball searches. The best detector for CEP, the AFS, was designed for high-p^ 
physics and added a CEP programme only in the final 2 — 3 years of operation; the analysis was 
completed after the ISR was turned off. The spectra for X = tt^-k' ,K^ K~ ,pp and 7r+7r~7r"'"7r~ 
showed structures (especially tt+tt") but no unequivocal evidence for glueballs. The low mass scalar 
sector is still unsettled, and the whole M^^r region from 400 — 1700 MeV/c^ may even be a very 
broad scalar glueball "cut" by the /o(980) and /o(1500). 

5 Fixed Target experiments 

After the ISR was closed in 1983, CEP of low multiplicity hadron states, for both hadron spectroscopy 
and for studies of production mechanisms, became an active field in fixed target experiments in the 
1980s and 1990s. We shall restrict ourselves to a few brief comments since our focus is on high 
energy colliding beam experiments. 

Thanks to high luminosities, dedicated experiments in multiparticle spectrometers, such as 
the Omega facility at the CERN SPS experiments (WA76, WA91 and WA102) and experiment E690 
at the Fermilab Tevatron, were highly productive. In particular, the experiments using the Omega 
spectrometer [114, 115] confirmed the earlier observation by the SFM collaboration (at the ISR) 
[116] of a correlation between the properties of the central system and the directions of the outgoing 
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protons. Close & Kirk suggested that glueballs and qq mesons may have different dependencies on 
the difference in the transverse momenta of the outgoing protons [117, 118] but this idea remains 
unproven. Moreover, s/s was only 29 GeV and so RIR and JRJP exchanges were certainly not 
negligible. Ref. [115] contains many relevant references to CEP at the Omega facility, with a large 
number of results that are of interest for meson spectroscopy. 

Experiment NA22 at the European Hybrid Spectrometer (EHS) also extracted DIPE cross 
sections in 7r+p and K^p collisions, although with pbeam = 250 GeV/c {y/s = 21.7 GeV) a large 
non-DFE background had to be subtracted [119]. They found an inclusive (but mostly tt+tt^) 
DIPE cross section of (39 ± 5 ± 8)/tib for K~^p and (24 ± 6 ± 3)y^b for Tr+p. Such a difference is 
qualitatively expected in Regge theory since ariK^p) > (jT(7r+p). 

At Pbeam = 800 GeV/c, Fermilab Experiment E690 studied the X = KgKg system, showing 
in particular a large /o(1500) peak [120]. 



In 1981 the first pp collider, the SppS, came into operation at CERN. It provided an increase in ^/s 
from 63 GeV at the ISR up to 630 GeV, and even 900 GeV for some low luminosity "ramping runs" . 
This step-up in energy enabled the discoveries of the W and Z bosons, and observations of dramatic 
high Et jets. Diffractive physics took a back seat during this period, although experiment UA8, by 
adding forward Roman pots on both beam pipes to the central UA2 experiment, observed diffrac- 
tivcly produced dijets [121] with Et up to ^ 13 GeV. Ingclman and Schlein had earlier proposed 
to describe the pomeron in hard processes as having quark and gluon constituents, analogous to a 
hadron [122]. In this picture, the parton densities in the pomeron can be measured in high mass 
diffraction. The jet kinematics gives the scattering parton momenta and thus rather direct informa- 
tion on the parton distributions in the IP. This approach led directly to the subsequent very precise 
measurement of diffractive parton densities at the HERA collider (see for example Refs. [8, 9]). 

Experiment UA8 measured p/p tracks after (low-/3) quadrupole magnets, and they studied 
p+X+p with X = "anything" , using the UA2 central detector [123]. They selected 107 events in the 
D IP E-dominated region with xp > 0.95 (i.e. gaps Ar] ^ 3 on both sides). Most of the events had 
Mx < 8 GeV/c^ but no specific exclusive states X could be resolved, as central particle momenta 
were not measured (UA2 had no magnetic field) and the missing mass resolution was much worse 
than the mass derived using the calorimeter alone (which was ctm ^ 1 GeV/c^). UA8 extracted the 
DIPE cross section, ajpp, assuming the validity of Eq. (2.9) with a flux (see Eq. (2.8)) given by 



6 CERN Proton-ant iproton Collider (SppS) 




(6.41) 
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and with parameters extracted from data on elastic scattering and single diffraction: 



a 



0.71 GeV^ 



b 



1.08 GeV 



-2 



K 



0.74 GeV 



-2 



a]p{t) 



1.035 + 0.165f + 0.059^2 . 



(6.42) 



UA8 inferred that app « 3 mb for 2 < Mx ^ 8 GeV/c^, decreasing to 0.5 - 1.0 mb for 
10 < Mx ^ 25 GeV/c^. This is significantly larger than expected assuming the simple factor- 
ization formula presented in Eq. (4.37) at the same total energy. The UA8 paper also included a 
prediction for daDFE/dMx at the Tevatron and LHC assuming a constant app = 1 mb. They 
predicted about 0.4 /ib/GeV/c^ at Mx = 100 (200) GeV/c^ at these colliders respectively. 

The largest experiment at the SppS, UAl, was designed to discover the W and Z, but it also 
performed a CEP study at y/s = 630 GeV using forward rapidity gaps [124]. They required no energy 
in forward calorimeters with 3 < jr^j < 6 and some energy in the central calorimeters (5° < 6 < 175°). 
A localized energy deposition of Et > 1-5 GeV and/or a "jet" with Et > 3 GeV was required 
to trigger, which excluded the possibility of detecting low mass exclusive states. The data were 
compared to a non-diffractive sample, containing at least one charged particle with 1.5 < \ri\ < 5.5 
in each forward direction. The mass Mx of the central hadrons had a mean (A-Ix) ^ 36 GcV/c^ and 
extended a little above 63 GeV/c^ (above which at least one scattered proton must have xp < 0.90). 
At low Mx, the mean charged particle multiplicity {rich) was measured to be similar to that in pp or 
e+e~ collisions at ^/s = Mx, but it was found to rise faster and, by Mx = 70 GeV/c^, it was a factor 
of two higher. This effect has not been quantitatively explained, and hopefully it will be checked at 
the Tevatron or LHC. Recall, however, that the final state in e~^e~ collisions producing dijets is qq 
and in IPIP collisions it is predominantly gg, so one might expect {richjiPiP > {nch)e+e- ^ ^ result 
of the difference between quark and gluons jets, together with pomeron "remnants". With the UAl 
jet algorithm it was found that 5% of the events have at least one jet with Et > 10 GcV, but this 
fraction would probably be lower if the trigger had not required an energy deposition of Et > 1.5 
GeV. This class of central dijet events is expected be a dominant background to searches such as 
p + H + p ^ p + JJ + p ai the LHC. 

7 Tevatron 
7.1 Introduction 

Until the LHC began operation in December 2009, the Tevatron at Fermilab was the highest energy 
hadron collider, making proton-antiproton collisions at y/s = 1.96 TeV. In Run I (1992-1996) it 

provided collisions at = 546 GeV, 630 GeV (to equal the CERN SppS collider) and 1800 GeV. 
The total delivered luminosity per experiment (CDF and D0) was 110 pb~^, compared to about 
7000 pb-i so far in Run II (2004-2009). 

In the first (Run I) phase of diffractive studies, installation of Roman pots by CDF, with 
silicon and drift chamber track detectors along both beam pipes, allowed measurements of elastically 
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and diffractively scattered protons and antiprotons. Elastic scattering was measured at ^/s = 546 
(1800) GeV for —0.29 < t < —0.025 GeV^, where it was well described by an exponential slope 
with b — 15.28 ± 0.58 (16.98 ± 0.25) GeV~^, and the integrated elastic cross sections were cxei = 
12.9 ± 0.3 mb and 19.7 ± 0.9 mb, at the two energies [125]. By measuring the rates of elastic 
scattering and inelastic collisions simultaneously, CDF exploited the Optical Theorem to measure 
the total cross section, aT, independent of the luminosity [126]. They found (Jt = 61.26 ± 0.93 
(80.03 ± 2.24) mb, with (Tci/o-tot rising from 0.210 ± 0.002 to 0.246 ± 0.004 over this y/s range. 
Two other experiments, E710 [127] and E811 [128] found significantly (~ 2.5c7) lower values for 
both C7ei and crtot (but with the same ratio) at ^/s = 1800 GeV. These results, together with LHC 
expectations, are discussed in [40] and [129]. CDF also measured the inelastic proton spectra for 
xp <^ 0.88 [130, 131], where they observed the expected dominance of the diffractive peak above 
Xi? ~ 0.95. This corresponds to diffractive masses Mx = 400 GeV/c^; indeed later Tevatron studies 
succeeded in observing W, Z and jets with Et up to 80 GeV in single diffraction. Unfortunately the 
forward proton detectors were removed without having measured CEP with both protons detected. 

The second stage of diffractive physics at the Tevatron took place in Run II, initially without 

any forward p/p measiircmcnts. CDF and D0 both stiidied hard single diffraction with forward 
rapidity gaps, after which CDF installed a new set of Roman pot detectors 50 m downstream to 
measure diffractive antiprotons. There were no detectors for forward protons, so p + X + p could 
not be studied with both the p and p tagged, however sets of scintillation counters with 9 < 3° were 
added along the beam pipes. These "beam shower counters" (BSC) acted as rapidity gap detectors, 
using showers produced in the beam pipe by secondaries with 5.4 < \r]\ < 7.4. When combined with 
"miniplug calorimeters" covering 3.6 < \r]\ < 5.2, gaps of Ar] > 3.8 could be required on each side, 
with or without a measured antiproton. Later D0 also installed Roman pots with tracking on both 
sides, but has not yet presented any CEP results. 

The single diffractive production of dijets [132, 133], W [134, 135] and Z [136] (preliminary) 
has been observed and the fraction of W, Z or dijets that are classed as diffractive, i.e. either with a 
high-XiT antiproton or a large forward gap, is around 1%. This is to be contrasted with ep collisions 
at HERA, where HI and ZEUS found that the corresponding fraction of diffractive to non-diffractive 
events is 10%. This is in accord with the rough expectation that the gap survival probability S"^ 
is lower by about an order of magnitude in hadron-hadron collisions than in 'y*p collisions. 

CDF also studied the central production of dijet events containing both a leading p with 
0.035 < < 0.095 and a rapidity gap in the proton direction, corresponding to 0.01 < < 
0.03 [137]. Although the proton is not observed, its ^ can be computed using 



where the sum is over all the observed particles in the detector. This constitutes an observation of 
DIPE dijet production, and a comparison to the theoretical predictions has been carried out using 
pomeron parton densities extracted from the HERA data in Rcfs. [138, 139]. Again a gap survival 
factor of around 10% was needed to fit the data. Importantly, it was noted that the rate of DIPE 
(two-gap) dijet events is approximately 10% of the rate of single diffractive (one-gap) dijet events , 
i.e. the price to pay for gap survival should not be applied twice. This was also observed in a CDF 
study [140] of inclusive DIPE, with a measured antiproton with 0.035 < < 0.095 and a rapidity 
gap on the proton side. 
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(7.43) 
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7.2 Central exclusive production 



The study of CEP at the Tevatron is a recent activity, largely motivated by the possibility of studying 
exclusive p+p p + H +p and related processes at the LHC. In 2001 there was a very large spread 
in predictions for this exclusive cross section, some of which even suggested that a Higgs search may 
be possible at the Tevatron, and a Letter of Intent to perform such a search by measuring both final 
state protons was submitted to the Fermilab Program Advisory Committee [141] . It did not proceed 
to the proposal stage however, partly because of the; large cross sc;c;tion uncertainty. Although the 
exclusive cross section is much smaller than the inclusive one, if the protons are measured precisely 
it was noted that the "missing mass" method can be used to give a very good measurement of 
the central state, independent of its nature [1]''. Although it is now understood that the Higgs is 
not likely to be produced this way at the Tevatron, the measurement of other central systems, in 
particular those with the same quantum numbers as the Higgs (such as XcO cind Xbo) constrains the 
predictions for a{p + H + p). The diagrams for XcO and Xbo production are similar to those for 
exclusive Higgs production, but with c- and 6-quark loops instead of a top-quark loop. Of course the 
relevant scale is much lower for these lighter states and this renders the applicability of perturbative 
QCD more uncertain. The cleanest related process is p + j^ + p (through IP IP — )• 77) but the cross 
section is small for St (7) > 5 GeV, where measurements were made. Exclusive dijet production has 
a much higher cross section but one must then compromise on not having such a clean final state. 




Figure 9: Predictions for CEP di-photon production at the Tevatron and the LHC. Figure from 
Ref. [142]. 



^Rcf. [1] also pointed out that a measurement of the difference in the proton arrival times at the forward detectors 
could be used to greatly reduce the impact of pile-up background. 
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Nevertheless, CDF has succeeded in measuring exclusive 77, Xc and dijct production, and we 
discuss each of these in turn. We also discuss exclusive 77 e+e^,/i+/i^ and 7JP J/ip,ijj'{2S), 
which was recently seen for the first time in hadron-hadron collisions. Apart from the dijet study, the 
forward p could not be detected, as the Roman pots (only on the p side) did not have acceptance for 
low Mx- The analyses instead relied on finding events that contain just the state X in conjunction 
with an otherwise empty (i.e. consistent with noise levels) detector. The miniplug calorimeters and 
BSC counters, which have coverage out to \ri\ = 7.4, were crucial for this exclusivity requirement 
(the BSC was used as a veto in the trigger) . 

Firstly we shall discuss the exclusive 77 search [143], which was combined with a CEP 
e+e~ search [144], as the trigger (and most of the analysis) is identical. Only in the final step was 
the central tracking used to separate 16 e+e~ events from three with electromagnetic showers with 
Et > 5 GeV and no tracks. In all cases the showers had tt and small, and the e+e" 

events agreed with the precise QED expectation, providing a good control for the 77 candidates. 
The gap survival probability is not an issue for the QED events; the impact parameter is large and 
5^ 1. Also, the balance in Et and Acf) should make it possible to find QED events (especially 
fi'^lj,~) in the presence of pile- up; this is now being studied in CDF. Background, e.g. from tt'^tt^ 
in the 77 candidate sample, could not be quantitively assessed a priori, but two of the events had 
narrow single showers on each side and were very unlikely to be background. The prediction using 

the Durham model [142] is shown in Fig. 9. The prediction of 86^3 fb for Exil) > 5 GeV and 
[77(7) I < 1 would give 0.8tJ;5 events, and the CDF data are in good agreement with this. The two 
events correspond to ^ lO"""^^ x (Xinei, showing that it is possible to find even very rare exclusive 
events. More CDF data has been taken with a lower threshold Et{j) > 2.5 GeV, and there are 
plans to search for exclusive 77 events at the LHC. 




Figure 10: Feynman diagrams for processes contributing to the exclusive di-lepton signal, (a) 
77 ^ 1+1-, (b) 7iP ^ J/V', V(25), ZO, and (c) FP ^ XcO- 

In addition to the exclusive 77 search, CDF also studied the production of lepton pairs 
(e+e~, /i+/i~), either in association with no other particles or with one additional photon. Such 

exclusive leptons may be produced through several mechanisms, as shown in Fig. 10. We begin by 
discussing 11^ fjT production at low M{n~^iJ,~), for which CDF used an exclusive di-muon trigger 
in the mass range M{ij,'^i2~) G [3.0,4.0] GeV/c^ and [jy^j < 0.6. The mass range was limited to 
M{fj,'^ljr) > 3 GeV/c^, as below pt = 1-5 GeV/c muons range out in the calorimeters. 

Photoproduction of a vector meson, shown in Fig. 10(b), is one source of exclusive muon 
pairs in both ep and pp{p) collisions. The predictions for the pp process are closely related to those 
for the ep process, bearing in mind the differing soft survival factors, S'^{pp) < S'^{ep). CDF recently 
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observed the production of J/^ and ip(2S), with subsequent decay to fi~^ ^~ , via this mechanism, 
the first such observation of vector meson photoproduction in hadron-hadron coUisions [145]. Simply 
requiring rapidity gaps in conjunction with two central muons results in an extremely clean spectrum, 
and any non-exclusive background (mostly undetected proton dissociation) is at most a few percent, 
see Fig. 11. The cross sections agree with predictions [146, 147, 148]. Although these vector mesons 
cannot be produced in either 77- or IP IP-coWisions, they could be produced in odderon-pomeron 
interactions. While the odderon is required in QCD [149], the couplings are expected to be suppressed 
compared to pomeron couplings and there is at present no direct evidence for it. 
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Figure 11: Di-muon mass distribution of exclusive fJ-'^fJ- events in CDF, with no EM shower, 
(histogram) together with a fit to two Gaussians for J/ip and ip{2S), and a QED continuum. Inset: 
Data above the J/ip and excluding 3.65 < M{fj,+fx-) < 3.75 GeV/c^ (V'(2S')) with a fit to the QED 
spectrum. Figure from Rcf. [145]. 



Another source of exclusive muon pairs, this time in association with a photon, is central 
exclusive Xc production: p+p^p + Xc+P, with subsequent decay: Xc — >■ J/''P + 7 ^ /^^M^ + 7j 
see Fig. 10(c). This is a process much more closely related, at least theoretically, to CEP of a 
Higgs boson. Measurement of this at the Tevatron provides, together with exclusive dijets which we 
shall discuss shortly, a good test of the theory. CEP of 77 and Xb are both superior theoretically, 
but unfortunately their cross sections are smaller. CDF has observed p + {xc J/ijj + j) + p 
events, with dcr/d?/|j^=o(Xc) = 76 ± 10 ± 10 nb, assuming the Xco dominates [145] (CDF could not 
resolve different Xc states). The theoretical predictions of the XcO cross section are in agreement 
with the measured value. Ref. [150] predicted d(T/d2/|j,=o(Xco) = 130 nb, however since that paper 
the PDG value of the XcO width was lowered by 45% [107], correcting their prediction to 90 nb 
and correcting Yuan's prediction [151] of 160 nb to 110 nb. Bzdak's prediction is 45 nb [152]. It 
should be noted however, that the assumed dominance of the XcO has been questioned. In Ref. [29] 
it was suggested that higher spin Xc states may contribute to the cross section, with their smaller 
production rates being compensated by larger branching ratios. Indeed, a recent re-analysis has 
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been performed by the Durham group and collaborators, including the effects of non-factorising, 
"enhanced", contributions to the gap survival probability [153, 73]. Their studies suggest that the 
XcO production cross section is somewhat lower than previous determinations, but that the reduction 
to the total rate is compensated by the production of Xci and Xc2 states. They find for the total 
production of all three states, da / dy\y=o{Xc) = 65 nb and the calculation has been implemented in 
the Monte Carlo SuperCHIC [154]. The agreement between theory and experiment is encouraging 
and indicates that the predictions for CEP Higgs production at the LHC are unlikely to be far ofi^ 
the mark. 

Exclusive muon [145] (and electron [144]) pairs are also produced via the QED 77 fusion 

process, depicted in Fig. 10(a). This process generates a continuum in the M(/z"'"/U~) distribution 
over the which the J/ip and ijj{2S) peaks are visible, as shown in Fig. 11. The continuum agrees in 
both shape and magnitude with the QED predictions. Also, since the exclusive p + l^l" +p cross 
section is very well known (only Coulomb elastic scattering is as well known), it is a good candidate 
for measuring the absolute luminosity of the LHC, and thus for calibrating high-rate luminosity 
monitors. The limit to the precision (^ 1% — 2%) with which this can be done will come from 
the knowledge of acceptances and efficiencies (favouring over e^er). One must do this in 

the presence of some pile-up, as the probability of vetoing on pile-up will depend on the unknown 
inelastic cross section. Another difficulty is that the protons may dissociate, e.g. p p* ^ Ntttt. 
The CDF BSCs, situated along both beam pipes, were efficient at detecting such dissociations, which 
create forward showers in the pipes and surrounding material. CMS and ATLAS do not yet have 
such counters, but there is a proposal to add them [155] . 

At high mass, lepton pairs may be produced in exclusive Z photoproduction, jlP — )• Z. 

A nearly real photon is radiated by one of the incoming beam particles, fluctuates into a qq pair, 
which then scatters by IP exchange on the other proton, see Fig. 10(b). CDF performed a search 
for this reaction [156], however, the cross section is predicted to be only 0.21 fb [157] or 0.3 fb [147] 
at the Tevatron, far below the CDF sensitivity. In contrast, as we shall discuss in section 8.5, the 
cross section at the LHC is expected to be about 6 fb for ly(2')] < 2 and the Z decaying to e+e" or 
fj.'^ fj~ [147]. It would be very interesting to measure exclusive Z production, as BSM theories with 
additional particles in the "flP — >■ Z vertex loop change the cross section. For example, in White's 
theory of the pomeron there should be a "large" (but not quantified) increase, due to colour-sextet 
quark loops [158]. Exclusive 77 pairs have been seen with pair masses up to 75 GeV/c^ and 

no observed Z candidates [156, 159]. 

Unlike exclusive lepton production, exclusive dijet production, X = jj, has the advantage 
of a relatively high cross section but, unlike AT = 77 or x^i it is subject to the usual uncertainties 
associated with defining jets. CDF measured it from a sample of DIPE dijet events, with a detected 
p together with a large gap on the p side [83]. They reconstructed the mass of the dijet system, 
Mjj, and the mass of all the particles in the central detectors, Mx- Exclusive dijets would have the 
ratio Rjj = Mjj/Mx, shown in Fig. 12, close to 1. Above Rjj ~ 0.6 the data show a significant 
excess over that expected from pomwig, which is a Monte Carlo generator for diffractive processes, 
including an implementation of the lugelman-Schlein model of DIPE [160]. pomwig is tuned to 
CDF and HI data and does not include an exclusive component. After combination with a full CDF 
detector simulation, the excess is well described by the predictions for exclusive dijet production 
given by the ExHuME Monte Carlo generator [161], which implements the Durham model. Making 
such fits for several values of the minimum jet transverse energy, E^™, CDF derived the exclusive 
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Figure 12: Dijet mass fraction for DPE data (points) and best fit to the sum of two components: 
POMWIG DPE + SD and ND background, and ExHuME (exclusive dijets). Figure from Ref. [83]. 



dijet cross section shown in Fig. 13. The data are far below the prediction of the dpemc Monte Carlo 
generator [162], which implements the Saclay model discussed in Section 3.2, and although slightly 
below the ExHuME prediction the agreement is much better and certainly within the theoretical 

uncertainties* . 

Exclusive dijets are expected to come primarily from gg — > gg since quark pair production, 
99 QQj is suppressed due to the = rule, by a factor ^ m^/Mjj [29]. In particular (see for 
example [30]) 

-^^99^99) = -^-^ (7.44) 

and 

-^{99 ^ 55) = ^1 - (7.45) 

where these partonic cross sections are to be used in conjunction with Eq. (3.22). Light quark jets 
cannot easily be experimentally distinguished from gluon jets, but some c-jets and &-jets can be 
tagged by secondary vertices or displaced tracks. CDF found the that fraction of all exclusive dijets 
that arc from c- or &-quark jets is suppressed at Rjj > 0.5 compared to lower Rjj values. This lends 
further support to the claim that the excess is due to CEP dijet production. Note that qq dijet 

suppression will be very important in reducing QCD backgrounds in exclusive H ^ bb searches at 
the LHC, as we discuss in Section 8.1. 

Tevatron luminosities are now (in 2009) too high, even at the end of a store, to study events 
with no pile-up. A search is being made in CDF [159] for exclusive T and high mass 



*We note that better agreement between the Durham approach and the CDF dijet data was found in the analysis 
of Ref. [163]. 
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Figure 13: Exclusive dijet cross section vs. the minimum Et of the two jets, compared to the DPEMC 
and ExHuME predictions. Figure from Ref. [83]. 

di-muons, with pilc-up, based on their distinctive kinematics (and no associated tracks on the /i^/i^ 
vertex) . Unfortunately Xfc T + 7 is not hkely to be visible since the cross section is expected to be 
much smaller than that of the Xc, and due to the difficulty in selecting the correct photon in events 
that are contaminated by pile-up. 

8 The Large Hadron Collider 

The idea to install detectors far from the interaction point at CMS and/or ATLAS with the capacity 
to detect protons scattered through small angles has gained a great deal of attention in recent years, 
and the report presented in Ref. [2] constitutes a significant milestone on the road to CEP physics 
at the LHC. Measurement of both protons by detectors located 420 m from the interaction point at 
ATLAS or CMS has the virtue that central systems with masses up to 200 GeV/c^ can be measured 
with an event-by-event precision of a{AIx) ~ 2 — 3 GeV/c^. Adding detectors at 220-240 m extends 
the acceptance to much higher central system masses, with the limiting factor becoming the rate of 
the production process. Fig. 14 shows the forward detector acceptance as a function of central system 
mass for the case where both protons are detected at 420 m and when one is detected at 420 m and 
one at 220 m. The clean environment of CEP, combined with four-momentum constraints, makes 
for reduced backgrounds (even in the presence of significant amounts of pile- up) , often aided by the 
fact that the centrally produced system is predominantly in a = 0, C-even, P-even state [29, 164]. 
Of course having such a spin-parity filter also provides an excellent handle on the nature of any 
new physics. For little extra cost, forward detectors promise to significantly enhance the physics 
potential of the LHC. 

That said, CEP is not without its challenges. As we have seen, the theory is difficult, 
triggering can be tricky, signal rates for new physics are often low, high precision tracking and 
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Figure 14: Forward proton detector acceptances, as a function of the Higgs mass (mass of the 
central system), in the case of two protons detected at 420 m (left) and in the case that one proton 
is detected at 420 m and one at 220 m (right). The different curves correspond to differing positions 
of the active edge of the tracking detector relative to the beam. Figure from Ref . [2] . 

timing detectors need building and installing and pile-up needs to be dealt with. As the last section 
showed, recent data on CEP from CDF at the Tevatron gives confidence in the theoretical modelling, 
and extensive studies have demonstrated that the other challenges can be met, e.g. see Ref. [2]. Both 
in ATLAS and CMS, detectors at 420 m are situated too far from the central detector to be included 
in the level 1 trigger, although this should be possible in a future upgrade of the data acquisition 
system. The data from the proton detectors are, of course, available for integration in a higher level 
trigger. Note that, unlike the 420 m detectors, those at 220-240 m could be included in the level 1 
trigger. 

8.1 Higgs: SM and MSSM 

The possibility of observing central exclusive Higgs boson production, p + p—>-p + H + p, is largely 
responsible for the current interest in exclusive processes at hadron colliders. 

There are several important features of the CEP process which may allow one to extract, 
potentially unique, information on the structure of the Higgs sector. Firstly, precision measurements 
of the scattered proton momenta, in very forward detectors, offer the possibility of observing the 
Higgs as a bump in the missing mass distribution [2]. This is possible with a resolution a{M) ~ 
2 GeV/c^ per event, regardless of the final state (provided one can trigger on the central system). 
The precise calibration of the forward detectors necessary to achieve such a measurement can be 
obtained from data, using another exclusive (QED) process: p + p ^ p + iJi^ i-r +p. The muons can 
be measured very well in the central detector, and then both proton momenta are known with an 
uncertainty limited by the incoming beam momentum spread, Ap/p ^ lO^**, i.e. Ap ~ 700 MeV. 
To get sufficient rate for calibration on a reasonable time scale (e.g. one day), one must use muon 
pairs with low mass, M^+^- > 10 GeV/c^. This is too small to have any acceptance for detecting 
both protons but by selecting same-side-dimuon events with \ri^± \ > 2, one proton can have large 
enough momentum loss to be detected. 

In addition to obtaining the mass and width (provided it exceeds a few GeV) of any observed 

Higgs, one can also probe its spin and CP properties, due to the J^*^ — 0++ selection rule. For 
non-zero angle scattering, the production of higher spin and odd parity states is allowed, though 
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they are strongly suppressed for high enough central masses by the smallness of the final state proton 
transverse momenta [164]. Given enough statistics one can distinguish the J = and higher spin 
states, by the azimuthal correlations of the protons. Recall that at the ISR (Section 4), one could 
separate X = tt+tt" data into S-wave (J = 0), P-wave {J = 1), and D-wave (J = 2) states, mass-bin 
by mass-bin, from the moments of the angular distributions of the pious, even finding a p peak that 
was too small to see in the raw Mj^+^- spectrum. This could be done at the ISR thanks to a high 
level of statistics that is not likely to be available in p+ p ^ p + bb+ p data at the LHC, but the 
principle is the same. 

The spin-selection rule also has the advantage that it leads to a leading-order suppression of 
the bb background in the CEP oi H ^ bb (the suppression goes like ^ nig/Mjj, see Eq.(7.45)). The 

production of exclusive dijcts from light qq is likewise negligible, but CEP of gluon pairs receives no 
such suppression. Exclusive dijets, X = J J, are thus mostly gg jets, with a modest admixture of 

bb jets. The gg dijets can be discriminated against by heavy-flavour tagging, using a combination 
of displaced vertices and leptons in the jets from both c- and &-quark decays. In the exclusive dijet 

case one can also profit from the fact that we can have JJ = bb or gg, but pairs with net fiavour 
(e.g. bg or be) are forbidden. Gluon splitting {g — )■ bb) can be reduced by requiring both jets to be 
tagged. 

Fig. 15 shows how the cross section for producing a CEP SM Higgs varies with Higgs mass 
for different gluon distribution functions. The cross section is small and leads to low production 

rates. In particular, the bb channel is very challenging, due to the aforementioned backgrounds 
and the additional background coming from pile up interactions at high luminosity, which we shall 
discuss in more detail later. Triggering in this case would certainly benefit from having 220-240 m 
detectors in place, but even then one relies on optimistic scenarios for the production cross section, 
detector acceptance and trigger efficiency. 

In contrast, the observation of a SM Higgs in the WW'-*'^ channel is more promising. An 
initial study of the fully leptonic {l~^l^^) and semi-leptonic {l^frrJJ) decay modes, with basic 
experimental cuts, was performed in [165]. The backgrounds considered in this study were 77 fusion 
and exclusive qq production with Ty-strahlung. The contribution of gg — ?• W^W~ through a quark 
loop was also considered, but found to be negligible. In order to bring these backgrounds under 
control a range of cuts were employed. The rapidities of the W decay products, in the signal case, 
are mostly central. In contrast, the backgrounds favour a more forward distribution, hence a cut of 
\irj\ < 2.5 was imposed. The 77 process is dominated by the region in which the outgoing protons 
carry very small transverse momenta, and a cut was therefore made to remove this region. Finally, 
the dijet invariant mass was constrained to lie within a window around the W mass, thus reducing 
the background from qq events not due to W decay. After imposition of these cuts one could expect, 
for Higgs masses in the range 140 GeV/c^ < mn < 200 GeV/c^, ~ 2 — 3 events in 300 fb~^ of 
data, with no appreciable background. One of these would be in the gold-plated fully leptonic 
channel, producing the striking signature of an event containing two leptons with no other tracks 

on the vertex and a large missing Et- A subsequent analysis of the Wiy'^*^ channel with a fast 
detector simulation, taking into account pile-up background and using a more refined triggering 
strategy, found that the signal should be observable, with signal to background ratio of 1 or better, 
for niH Z 120 GeV/c^ and 300 ib'^ of data [2]. 

Unlike the SM case, the bb channel can be rather easier to explore in certain MSSM scenarios. 
The neutral part of the MSSM Higgs sector consists of one light and one heavy CP-even scalar, the 
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Figure 15: CEP cross section for Standard Model Higgs production. Figure from Ref. [2]. 

h and H respectively, and one CP-odd scalar, the A. For large tan /? and small niA, the hbb coupling 
is strongly enhanced over the SM case. Fig. 16 shows an example of the region of parameter space^ 
in which one could observe h ^ bb using CEP, for different amounts of integrated luminosity. Fig. 17 
shows the result of an in-depth analysis of one particular point in the niA ~ tan f3 plane (tan /3 = 40 
and ruA — 120 GeV/c^) [166]. The details of the two analyses can be found in Refs. [2, 166, 167, 168] 
and they are in broad agreement. 

A similar pair of plots for H ^ bb are shown in Fig 18. In addition to sensitivity in the low 
niA region, one can also hope to probe reasonably large values of tua (at large enough tan (3). If the 
MSSM parameters are in this region, an observation using CEP could be critical to understanding 
the Higgs sector at the LHC. This is because, in the large tua limit, referred to as the "decoupling 
limit", the light CP-even Higgs becomes SM-like. In contrast, the coupling of the H and A (which 
are approximately mass degenerate) to SM gauge bosons is strongly suppressed. As a result, the 
dominant decay of the H and A is to bb (r+r" typically contributes approximately 10%). This leads 
to the so called "LHC wedge region" [169, 170, 171], where the H and A may escape detection in 
inclusive searches. Not only that, but current inclusive techniques, proposed to probe the spin and 
CP properties of the Higgs sector, rely on either the ZZ decay mode at high mass [172] or vector 
boson fusion at low mass [173, 174, 175], both of which would be suppressed in this case (see also 
the discussion of this point in Ref. [168]). It should also be noted that this decoupling of higher mass 
Higgs states from the SM gauge bosons is a common feature of extended Higgs sectors, since such 
couplings are constrained by the electroweak precision fits. Furthermore, the CEP production of the 
A boson is strongly suppressed by the parity-even selection rule. This means that any measurement 
of the H mass and width using the forward proton detectors is unaffected by a nearby A. This is 
not the case in inclusive production. 

^In the M™^'" scenario with ^ = -1-200 GeV. 
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Figure 16: 5a discovery contours (upper plot) and contours of Scr statistical significance (lower plot) 
for the h ^ bb channel in CEP in the Ma - tan/3 plane of the MSSM within the Mh^max benchmark 
scenario (/i — +200 GeV) for different luminosity scenarios. The values of the mass of the light 
CP-even Higgs boson, Mh, are indicated by contour lines. The blue shaded region corresponds to 
the parameter region that is excluded by the LEP Higgs boson searches, while the purple shaded 
region is that excluded by searches at the Tevatron. Figure from Ref. [168]. 
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Figure 17: Upper: A typical mass fit for the 120 GeV/c^ MSSM h bb after 3 years of data 
taking at 10'^'* cm"'^ s~^ and after removing the overlap background contribution completely (e.g. 
using improved timing detectors). The significance is 5<t for these data. Lower: Significance of the 
measurement of a 120 GeV/c^ MSSM Higgs boson versus luminosity, for two different combinations of 
muon (MUG, MUlO) and fixed-jet-rate (FRT25, FRTIO) triggers. In each case the curves correspond 
to the baseline (lower) and improved-timing (upper) scenarios, i.e. the curves labelled "OLAP" 
include pile- up background. Figure from Ref. [2]. 
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Figure 18: 5a discovery contours (upper plot) and contours of Scr statistical significance (lower plot) 
for the H ^ hh channel in CEP in the Ma - tan/3 plane of the MSSM within the benchmark 
scenario = +200 GeV) for different luminosity scenarios. The values of the mass of the heavy 
CP-even Higgs boson, Mh, are indicated by contour lines. The blue shaded region corresponds to 
the parameter region that is excluded by the LEP Higgs boson searches, while the purple shaded 
region is that excluded by searches at the Tevatron. Figure from Ref. [168]. 
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Returning to Fig. 17, the curves in the lower plot correspond to diflfercnt trigger scenarios. 
They also indicate the influence of pile-up background, which is of concern to any CEP analysis. 
For brevity we denote a signal event with two detected protons and a central state all arising from 
the same collision as [pXp]. Pile- up backgrounds can originate when one or both protons and the 
central system originate from different pp collisions. If they originate from three different collisions 
(triple pile- up) wc write [p] [X] [p] and if they originate from two different collisions (double pile-up) we 
write [pX] [p] or [pp] [X] . In the latter case, since there is no acceptance for pure elastic scattering, the 
central detector will contain particles from both the collision that led to the two measured protons, 
as well as the particles from the second collision (which has no detected protons). The double and 
triple pile-up rates both grow strongly with luminosity and &i L k, 10^'* cm^^s^^ [p][X][p] typically 
dominates [2]. Fortunately there are many ways of reducing the pile- up background to manageable 
proportions. By measuring the four-momentum of the protons, one can determine the mean rapidity, 
Pt and invariant mass of the central system. Consistency between these values and those obtained 
using the central detector then provide a strong constraint. Furthermore, a very large reduction in 
the pile-up background can come from a precision measurement of the difference in time At between 
the arrival of the two protons in the forward detectors [2]. If they came from the same collision, 
displaced a distance z^p from the nominal vertex, then At ~ 2zpp/c and one can require agreement 
between the measurements of Zpp by the forward proton detectors and as determined from the central 
tracks. The resolution is (j{zpp) = 2.1 mm (3 mm/-\/2) for cr(At) = 10 ps. As the distribution of 
interactions has a spread around 60 mm (at la) this leads to a large reduction in pile-up background. 
Detectors with time resolution ^ 10 ps have been developed for this purpose, and there is hope for 
further improvement (the spread in arrival times due to different proton trajectories is only ~ 2 ps). 

In principle there are two other ways in which pile-up background can be reduced by fast 
timing. Firstly, additional calibrations could permit the actual collision time to be measured using 
the forward protons. The collisions in a bunch crossing have a time spread ^ 160 ps. If the 
individual collision times could also be measured using the detectors in the central region, a match 
with the collision time from the protons would give additional pile-up rejection. At present the 
central detectors do not have sufficient time resolution to make this practical. Secondly, one could 
place arrays of small (^ 1 cm^) fast-timing detectors in the forward regions, e.g. ^ < \rj\ < 5. 
Most inelastic collisions will produce particles in both regions, but central exclusive events will not. 
One can then veto events that contain extra particles originating from the collision vertex (which is 
defined by the measured protons). In the limit of the complete elimination of pile-up background, one 
could make a 5(7 discovery of an MSSM Higgs boson with 3 years of high luminosity (10^^ cm~^s~^) 
data taking (the mass peak is illustrated in the upper plot in Fig. 17). 

8.2 Higgs: NMSSM 

We now look at the capabilities for CEP of Higgs boson production in the Next-to-Minimal Super- 
symmetric Standard Model (NMSSM). More details can be found in [176]. Before discussing the 
specific details, let us emphasise that this scenario is particularly interesting since it deals with the 
possibility that a Higgs of mass 100 GeV/c^ decays predominantly to light particles (in the case 
we shall review it is four taus). Such a troublesome decay channel should be taken seriously and, as 
we shall see, CEP can provide an excellent means to study decays of this type. 

The NMSSM is an extension of the MSSM that solves the /^-problem, and also the little 
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hierarchy problem. It achieves this by adding a gauge-singlet superfield, S, to the MSSM such that 

the /i term is now dynamical in origin, arising when the scalar member of S aquires a vacuum 
expectation value (VEV). Since /x is no longer fundamental and therefore is no longer naturally 
of order the GUT scale (as would be the case if it were the only dimensionful parameter in the 
superpotential) , the /U-problem is solved. The little hierarchy problem is also solved because a 
lighter Higgs is allowed, thereby taking the pressure off the stop mass. mj. More specifically, the 
lightest scalar Higgs can decay predominantly to two pseudoscalar Higgses and the branching ratio 
to 6-quarks is correspondingly suppressed, thereby evading the 114 GeV/c^ bound from LEP (it 
drops to 86 GeV/c^). Having a lighter Higgs means that mj does not need to be so large, and that 
is preferred by the Z boson mass. 

The Higgs sector of the NMSSM extends that of the MSSM by adding an extra pseudoscalar 

Higgs, the a and an extra scalar Higgs; crucially 5 is a gauge singlet and hence h ^ aa can dominate 
with a light a (i.e. below the threshold for a — >■ bb). Freed of the heavy i, it is natural to have a light 
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Figure 19: CEP of an NMSSM Higgs. 



Higgs with a reducing branching ratio to &-quarks. Masses of the h in the range 85-105 GeV/c^ are 
most natural in this scenario and our attention will focus on ruh = 93 GeV /(? and nria =9.7 GeV/c^ 
with BR(/i aa) = 92% and BR(a tt) = 81% [177]. The hglitness of the pseudoscalar a 
means that the h decays predominantly to At. Should such a decay mode be dominant at the LHC, 
standard search strategies could fail and, as we shall see, CEP (as illustrated in Fig. 19) could provide 
the discovery channel. This "natural" scenario of the NMSSM has two additional bonus features 
that one might draw attention to: (1) A light Higgs is preferred by the precision electroweak data 
(recall that the best fit value is somewhat below 100 GeV/c^), and (2) A ^100 GeV/c^ Higgs with a 
reduced (~10%) branching ratio to 6-quarks naturally accommodates the existing 2.3 a LEP excess 
in e+e- Zbb [178, 179]. 

To detect the 4t decay of an NMSSM Higgs using CEP, we need first to trigger on the event. 
In this regard, we require that at least one of the rs decays to a sufficiently high px muon (or two 
decay to lower px muons). The detailed analysis is outlined in [176], here we shall just highlight 
the key features. Table 2 shows how the signal (CEP) and backgrounds are affected by the cuts 
imposed. The backgrounds are typical of CEP Higgs studies: DPE refers to the DIPE production 
of a pair of 6-quark jets (which is the dominant DIPE background after 6-tagging both jets) and it 
was simulated using the pomwig Monte Carlo [160]; OLAP refers to the background due to pile-up, 
in particular the three-fold coincidence of two single-diffractive events with an inclusive pp ^ X 
event. The QED background arises from 77 — > 4t and is easily removed. 



40 




M (GeV) 



Figure 20: A typical a mass measurement. Figure from Ref. [176 



The top line of the table is the cross section after requiring that there be at least one muon 
with pt > 6 GeV/c, which is the nominal minimum value to trigger at level 1^*^, and the condition 
that both protons be detected in the 420m detectors. There is also a loose cut on the invariant 
mass of the central system. Of the remaining cuts, we single out the "iVch = 4 or 6" cut for special 
mention. The charged track (Nch) cut is noteworthy because it can be implemented at the highest 
LHC luminosities: we cut on 4 - 6 charged tracks that point back to the vertex defined by the muon. 
Pile-up events add extra tracks (to both the signal and background) , but they do not often coincide 
with the primary vertex (e.g. within a 2.5 mm window) and do not spoil the effectiveness of this cut. 
The ability to make such hard cuts on charged tracks is of much wider utility than in this analysis. 
The 4-6 track event is then analysed in terms of its topology, exploiting the fact that the charged 
tracks originate from four taus in the signal, which themselves originate from two heavily boosted 
pseudoscalars. This analysis is purely track-based with almost no reliance on the calorimeter, which 
is more affected by pile-up. 

Accurate measurements of the proton momenta allow one to constrain the kinematics of the 
central system (in particular its invariant mass, px, and rapidity are known). We can also extract the 
masses of the h and a on an event-by-event basis. The mass of the h is straightforward (it is measured 
directly by the forward proton detectors) and a precision below 1 GeV/c^ can be obtained with just 
a few events. The measurement of the pseudoscalar mass is more interesting and potentially very 
important. The proton measurements fix p^, and Px.y ~ for the central system. In addition, the r 
pairs are highly boosted, which means they are coUinear with their parent pseudoscalars. Thus the 
four-momentum of each pseudoscalar is approximately proportional to the observed (track) four- 
momentum. The two unknown proportionality constants (i.e. the missing energy fractions) are 
overconstrained, since we have three equations from the proton detectors. The result is that we 
can solve for the pseudoscalar masses, with four measurements per event. Fig. 20 shows a typical 
distribution of a masses based on 180 fb~^ of data collected at 3 x lO'^'^ cm^-^s^^. 

Although we can expect only a small number of signal events, the background is under 

^''it will turn out that a higher cut of 10 GeV/c is preferred in the subsequent analysis. 
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CEP 


DPE 


OLAP 


QED 


Cut 


H 


bb 


99 


bb 


bb 


At 


2t 21 


pt^o, a, 6, M 


0.442 


25.14 


1.51x10^ 


1.29x10^ 


1.74x10*5 


0.014 


0.467 


TVch = 4 or 6 


0.226 


1.59 


28.84 


1.58x10^ 


1.44x10^ 


0.003 


0.056 


Qh = 


0.198 


0.207 


3.77 


18.69 


1.29x10^ 


5x10-^ 


0.010 


Topology 


0.143 


0.036 


0.432 


0.209 


1.84 




<0.001 


Py, isolation 


0.083 


0.001 


0.008 


0.003 


0.082 






Pt 


0.071 


5x10--* 


0.004 


0.002 


0.007 






nia > 2mr 


0.066 


2x10"^ 


0.001 


0.001 


0.005 







Table 2: The table of cross sections for the signal and backgrounds for NMSSM Higgs production. 
All cross sections are in femtobarns. The pile-up (OLAP) background is computed at a luminosity 
of 10^^ cm-2 s-^ Prom Ref. [176]. 

control. One should be able to collect 4 signal events, on a background of 0.1 events, with 3 years 
running at an instantaneous luminosity of 5 x 10^^ cm~^s~^. Even with these few events, one may 
extract the masses of both the a and the h, with a resolution below 2-3 GeV/c^. The statistical 
significance of any discovery can exceed 5cr in a few years of running at luminosities in excess of 
10^3 cm-2s-i. 

8.3 Other new physics 

In addition to the MSSM and NMSSM Higgs studies, the CEP of a number of other BSM central 
systems has been considered in the literature. 

For instance, in [180] a CP-violating version of the MSSM was studied. In such a model, 
the h, H and A may all mix to produce mass eigenstates. Hi, H^, of indefinite CP. For large 
tan/3, the couplings of the Hi to bb and r+r^ are enhanced, the states are naturally nearly mass 
degenerate (separated by a few GeV/c^) and strong mixing between the states leads to a complicated 

resonance structure. The CEP of the Hi, with decay to bb and t+t~ was studied and it was shown 
that the missing mass distribution could be used to perform a line-shape analysis of the resonance 
structure, which is probably impossible in inclusive production. Also in the CP-violating MSSM, but 
in a different parameter region, the LEP bounds on the lightest Higgs mass become much weaker. 
In [181, 182], CEP of such a light Higgs was studied, though in this case the bb channel appears not 
to be observable and the r+r" channel is marginal. 

Another extension of the Higgs sector studied in relation to CEP is a triplet Higgs model 
in which the Higgs sector consists of one complex Higgs doublet and two triplets (one real and one 
complex) [183]. In such a model, the neutral Higgs associated with the doublet, has an enhanced 
coupling to fermions^^. The observation, using forward proton detectors, in the H^ — )• bb channel 
was considered and it was found to be observable, with a significance ranging from ~ 4cr — 14cr over 

11 The becomes the Standard Model Higgs boson for vanishing doublet-triplet mixing. 
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much of the parameter space with 60 fb ^ of data. With such data, the mass could also be measured 
very accurately, with an expected error less than 2 GeV/c^ over a significant region of parameter 
space (and in some regions better than 0.3 GeV/c^). 

In some extensions of the SM the Higgs decays predominantly to unobservable particles. 
The CEP of such an "invisible Higgs" has been studied in [184]. However, while a detection in the 
CEP channel would be especially useful, since it would not rely on the detection of the Higgs decay 
products, it is not clear that pile-up backgrounds can be brought sufficiently under control. 

CEP Higgs production has also been studied in a model where the SM is complemented by 
a fourth generation of very heavy fermions [168]. The Higgs coupling to gluons in such a model is 
typically larger than in the SM. In the H ^ bb channel, it was found that a significance in excess 
of 3(7 could be obtained with 60 fb~^ of data over the range 114 GeV/c^ ^ mn ^ 145 GeV/c^ not 
excluded by LEP or the Tevatron^^. 

Away from the Higgs sector, the production of long-lived gluinos, expected in the "Split 
Supcrsymmetry" model [185, 186], was studied in Rcf. [187]. Such particles hadronisc to produce 
so-called i?-hadrons which look much like slow muons. The CEP process in this case is particularly 
clean, since cuts on the velocity of the i?-hadrons reduce the muon background to a negligible level. 
It was found that the measurements from the forward proton detectors could in this case give a 

measurement of the gluino mass to better than 1%, for masses below 350 GeV/c^ and with 300 fb~"^ 
of data. Although not a discovery channel, such a measurement would be highly complementary to 
any measurement in inclusive production, which suffers from large systematic uncertainties in this 
mass region [188]. 

Another model relevant to CEP is White's theory of the critical pomeron [158], which 
requires the existence of a pair of massive colour sextet quarks, Qg. Since Qs loops couple strongly 
to the pomeron and to weak bosons, processes such as PIP — > W^W~ and ZZ (but not WZ) should 
occur with a much enhanced cross section. The threshold for this process is above the Tevatron 
energy but below the LHC energy. The enhanced cross section may even be visible without forward 
proton tagging, but the latter will be needed to demonstrate that PIP interactions are responsible. 
Photoproduction of Z, 'yP — )• Z, should also be enhanced a factor > 10 over the SM, and should be 
detectable. 

8.4 Two photon collisions 

So far in this review our attention has focussed mainly on CEP via the strong interaction. At the 
LHC, high proton luminosities will ensure that the flux of bremsstrahlung photons off the colliding 
beams is sufficiently large for CEP via 77 — )■ X and jP — > X to be of interest. We focus in this 
section on the former. The elastic photon flux is shown in Fig. 21, from which the corresponding pp 
cross section, dupp, can be obtained: 



where the integral is over the 77 centre-of-mass energy, W^^. To put this flux in context, we note 
that it corresponds to a 77 luminosity that is 1% of the pp luminosity for the production of central 

^^Thc region mn ^ 220 GeV/c^ is also experimentally allowed, however the CEP rate for H bb is too low to be 
accessible in this region. 




(8.46) 
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Figure 21: Luminosity spectrum for photon-photon coUisions at the LHC in the range Q^j„ < < 
2 GeV^ (sohd blue hue) compared to the corresponding luminosity if the energy of each photon is 
restricted to the forward detector tagging range 20 GeV < < 900 GeV (dashed green curve). 
Figure from Ref. [189]. 

systems with mass greater than 23 GeV/c^, falling to 0.1% for masses above 225 GeV/c^. In addition, 
the photon virtualities are cut-off sharply by the proton electromagnetic form factor and are small, 
i.e. (Q^) « 0.01 GeV^; they can usually be assumed to be zero. 

Our discussion in this section follows closely the presentation in [2, 189]. In turn, we will 
consider CEP of: (i) Icpton pairs; (ii) weak vector boson pairs; (iii) supersymmetric (SUSY) pairs. 
The first of these (especially using muon pairs) is a good candidate for measuring the LHC absolute 
luminosity, the second is an accurate probe of the ^^WW coupling and the latter could be used to 
make accurate mass measurements of the SUSY particles from the forward protons. 



Processes 


a(fb) 


Generator 


77 ^ (p^ > 2 GeV/c, [t?^! < 3.1) 


72 500 


LPAIR [190] 




108.5 


mg/me [191, 192] 


F+P- [M = 100 GeV/c^) 


4.06 


// 


F^F- {M = 200 GeV/c2) 


0.40 


// 


S+S- (M 100 GeV/c2) 


0.68 


// 


S+S- (M = 200 GeV/c2) 


0.07 


// 


H ^bb (M = 120 GeV/c2) 


0.15 


mg/me [191, 192] 



Table 3: Production cross sections ioi pp ^ p + X + p (via 77 exchange) for various processes {F = 
fermion, S = scalar). M is the mass of the corresponding particle. Table from Ref. [189]. 

Using a modified version of MADGRAPH/madevent [191, 192], the Louvain group [189] has 
computed the production cross sections for the various processes of interest. These are tabulated in 
Table 3. Since the cross sections for pair production depend only on the charge, spin and mass of the 
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Figure 22: Cross sections for various 77 processes at the LHC as a function of the minimal 77 
centre-of-mass energy Wq. Figure from Ref. [189]. 

produced particles, the results are shown for singly-charged, colourless fermions and scalars of two 
different masses. The cross sections are also shown as a function of the minimal 77 centre-of-mass 
energy Wq in Fig. 22. Note the hierarchy of cross sections, which is driven by the spin of the relevant 
t-channel exchange. The cross sections are large enough to merit further study. The two-photon 
production of SUSY pairs was first studied in [193] and the exclusive two-photon production of a 
SM Higgs was studied in Refs. [194, 195]. The rate is significantly smaller than that arising from 
strong dynamics (see Fig. 15) but it could become particularly interesting in the case of an enhanced 
i?77 coupling. 

Lepton pairs 

Two photon collisions can produce exclusive pairs of any charged particles, but it will probably not 
be possible to observe 77 — >■ qq because of overwhelming strong interaction backgrounds {DPE). 
However, lepton pairs (e^e^,/i^/i^ and also, to some extent, r+r^) should have very small back- 
grounds, even in the presence of pile-up. One requires that there are no additional charged particles 
on the dilepton vertex, and that pt{1'^1~) is very small or, equivalently, that the leptons are highly 
back-to-back ([tt — A(/)| <^ 1.0 GeV/Af;+;- rad) and with similar px- The A0 cut is especially pow- 
erful as it is forgiving of poor momentum resolution at high mass. The cross sections, calculated by 
the LPAIR Monte Carlo [190] at y's = 14 TeV, integrated over M^^ > Mi^i_^{min) and for different 
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are well approximated by the purely empirical fit: 



(r{Mmin,Vmax) = (480 - 1.25 X Mmin) X 




pb. 



(8.47) 



This equation agrees with the predictions of LPAIR to within ±4% for M„j„ from 8 GeV/c^ 
to 100 GeV/c^ and out to rjmax = 3. This pure QED cross section is so well known that it may 
provide the most precise measurement of the luminosity of the LHC [195, 196, 197, 198, 199, 200, 201] 
(integrated over a period), and help calibrate the luminosity monitors, which give a prompt rate 
measurement. For example the cross section for M^+^- > 10 GeV/c^ and \r]\ < 2.0, where one should 
have excellent and very well-known acceptance, is 8.55 pb. In 500 pb~^, with an unprescaled trigger, 
one would therefore have ^ 4275 events, which translates into a ~ 1.5% statistical uncertainty. A 
knowledge of the acceptance and efficiencies at this level is thus desirable. Without detecting both 
protons (and no LHC experiments will be able to do that in normal low-^ (high luminosity) running 
at such low masses) a concern is that one or both protons may dissociate. This probability can be 
measured in shower counters along the beam pipes at low luminosity, when pile-up can be neglected. 

Although the QED calculation of the dilepton cross section is accurately known, one still 
needs to worry about the issue of gap survival S"^. We have already seen that this suppression is 
large for CEP Higgs production, however 77 collisions are much more peripheral and the gap survival 
factor is much larger. Nevertheless, it must be understood at the percent level if we arc to exploit 
CEP dilepton production as a means to measure the luminosity. Fortunately, studies indicate that 
the gap survival should be very close to 100% for low-mass dilepton production [201]. 

How does 77 — )■ r+T^ compare with exclusive p + p — >■ p + H + p — >■ p + t^t^ + p? 
Differentiating the above equation we get da/dM w 0.4 fb/GeV/c^ with |?7(r)| < 2.0 at 120 GeV/c^. 
The branching fraction for a SM Higgs (120) H r+r" is ~ 0.07 and so, if a{p + H+p) ~ 3 fb [30], 
the rates will be at least similar. For example, in a 4 GeV/c^ mass window (which is commensurate 
with the expected mass rcsohition of the forward detectors) the cross sections would be aromid 1.6 fb 
and 0.2 fb respectively. Even assuming no other backgrounds and 100% detection efficiency, this 
translates into no more than 2a evidence of H ^ t~^t~ with 100 fb~^ of data. It will therefore be 
difficult to see a SM Higgs in this decay channel at the LHC, but it may be feasible if the Higgs 
production rate is enhanced (e.g. see [168]). The 77 process is a continuum with well-known cross 
section on which H — > t+t^ would (if the production cross section is sufficiently large) be a narrow 
peak. By measuring the protons, one could exploit the fact that the background is peaked at smaller 
values of \t\ than the signal to improve the significance of any observation. 

As discussed earlier, exclusive events are also very important for the absolute calibra- 

tion of the forward proton spectrometers. Note that the two-photon exclusive production of e+e" 
pairs can also be studied at the LHC, though triggering of such events is more difficult, and the Acp 
and pT(e"'"e~) measurements are worsened by bremsstrahlung. 

Finally, we note that the prospect of studying so called "unparticles" using the exclusive 
77 — )• channel has been studied in [202]. 
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W and Z boson pairs 



As illustrated in Fig. 22, the rate for pair production of W bosons via 77 — >■ W^W~ is large (around 
100 fb) and this can be exploited to investigate the corresponding quartic coupling. The signal is 
very clean and offers the possibility to study this process out to Mww ^ 500 GeV. The cross section 
for > 10 GcV/c and < 2.5 bom p + p -^p+ijj W+W' -> ^i+^i-vD) +p is 0.76 fb [189]. 
The large reduction comes from the two W ^ fiv decays. The cross section is only slightly reduced 
after adding the requirement of at least one forward proton tag. The unique signature of W pairs in 
the fully leptonic final state, no additional tracks on the l^lr vertex, large lepton acoplanarity and 
large missing transverse momentum strongly reduces the backgrounds. At low LHC luminosities (i.e. 
■C 10'^^ cm~^s~^), pile-up is not a problem and the events can be recognised even without detecting 
the outgoing protons [189, 203]. However tagging the protons would give improved sensitivity and a 
M\vw measurement. At higher luminosities, proton tagging is needed in order to suppress pile-up 
backgrounds. The situation is similar for the production of Z boson pairs, assuming fully leptonic, 
or semi- leptonic decays. In the SM, 77 — >• ZZ is negligible (it is absent at tree level). Thus any 
observation would be a direct measurement of an anomalous 'jjZZ coupling. 

In CEP olW~^W~ and ZZ, the fully leptonic decays (e.g. WW -> e~ ii'^fr) have practically 
no background, but constitute only 4.4% and 0.45% respectively of all decays. However in CEP one 
may be able to use additional four-momentum constraints to perhaps use about 50% of the events, 
i.e. all decays except the (six-jet) fully hadronic. For example, in p+p — ^ p+WW+p — >■ p+lvJJ+p, 
take the four-momenta of all detected objects p,p,l,J,J and calculate the missing mass, which is 
m„ 0. Similarly if X = ZZ l'^l~uD the missing mass to p,p, is mz (provided the missing 
Z ^ vD is on-shell). CEP can therefore produce a mass peak in Z ^ vv. With similar 4-momentum 
constraints as much as 50% of all (on-shell) W^W~ and ZZ events can be used. 

Anomalous 'y^WW couplings were first studied, within the context of a future e+e~ collider 
in [204]. Under the assumptions of electromagnetic gauge invariance, C and P conservation and 
custodial symmetry (to keep the tree-level p = M^/(Mf cos^ 9w) parameter equal to unity) there 
are only two relevant couplings to lowest order in an expansion in the inverse of the scale of new 
physics. A, i.e. the lowest non-renormalizable terms occur at dimension 6 in the lagrangian [204]: 

2 2 
Le = -j^F^.F^^'-W+'^W- - ^^F^^F>'^{W+"W0 + W-"W+). (8.48) 

The anomalous couplings are denoted by ao,c and they are zero in the Standard Model. Sensitivity 
to these anomalous couplings has been investigated by the Louvain group [189]. They focussed 

upon 77 — W~^W~ — >■ l^l~vv and 77 — ZZ using the signature of two leptons (e or 

/x) within the acceptance cuts jyyj < 2.5 and px > 10 GeV/c. With just a few fb~^ of data, the 
limits are expected to be orders of magnitude better than the best limits established at LEP2 [205] 
and significantly better than those that would be obtained using the W77 final state at the LHC 
[206, 207]. Similar conclusions have been reached in the studies presented in [208]. 

A study of the anomalous triple gauge couplings has also be performed. However, in this 
case the expected sensitivities are not as impressive [209]. 
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Supersymmetric pairs 




7 





7 





Figure 23: Relevant Feynman diagrams for SUSY pair production with leptons in the final state: 
chargino pair production (left); slepton pair production (right). Figure from Ref. [189]. 

The LHC, of course, is geared up to discover and explore low-scale supersymmetric particle 
production. However, SUSY pair production in 77 collisions, with two tagged protons, offers the 
opportunity to make direct measurements of the sparticle masses [189, 210]. As we shall see, the 
cross sections are not large and sufficient rate is only possible if the sparticles are light. The two- 
photon production of supersymmetric particles has been investigated in [193, 211, 212]. For chargino 
or slepton pair production, the final states are simple, see Fig. 23. Selecting a final state composed of 
two charged leptons with large missing energy and lepton acoplanarity means that the backgrounds 
should be under control, and 77 — >■ W^W~ generates the only irreducible background. 

In Ref. [189], three benchmark points in the parameter space of the constrained minimal 
supersymmetric model (CMSSM) were studied^'^. Here we will present the results for point LMl, 
which predicts a 97 GeV/c^ neutrahno (the LSP), a of mass 118 GeV/c2 (? = e,At), a [J of mass 
184 GeV/c^ and a lightest chargino, Xi 1 of mass 180 GeV/c^. Signal and background samples were 
produced using a modified version of CALCHEP [214] and the following acceptance cuts were applied: 
two leptons with > 3 GeV/c or 10 GeV/c and \rj\ < 2.5. Standard high level trigger efficiencies 
are high for all these types of events. 

Measuring the outgoing protons allows one to measure the cross section as a function of M^^, 
and expected event rates with 100 fb~^ of data are presented in Fig. 24. The two edges, at twice 
the mass of the and at twice the mass of the are visible. However, the event rates are low and 
mass measurements with a precision of a few GeV/c^ will require very high statistics. HECTOR [215] 
simulations of forward protons from slepton events consistent with the LMl benchmark point indicate 
that the 220 m detectors will have both protons tagged for only 30% of events. Addition of detectors 
at 420 m increases that to 90% of events. 

Another attractive feature of two-photon production of SUSY pairs with tagged protons 
is that both the missing energy f and fjx are measured, as well as the missing mass (here called 
W^miss)- For the SUSY signal W^\ss starts at 2 x mi5p[210], as illustrated in Fig. 25, while the 
background starts at Wmiss = 

^•^The parameter space is constrained to be in agreement with the latest cosmological data [213]. 
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Figure 24: The invariant mass spectrum for SUSY pair production in the LMl scenario and assuming 
100 fb~^ of data. The WW background has been down-scaled by the quoted factor. See text for 
event selection criteria. Figure from Rcf. [189]. 
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Figure 25: The missing mass distribution for SUSY pair production in the LMl scenario and assum- 
ing 100 fb^^ of data. Event selection assumes two isolated opposite charge muons with pr > 7 GeV/c 
or two electrons with pT > 10 GeV/c, < 2.5 and two tagged protons. Figure from Ref. [210]. 
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8.5 7P collisions 



In addition to DIPE and two-photon collisions, it is also possible to study exclusive photoproduction 
(7JP) at the LHC. A quasi- real photon is radiated off one of the protons and fluctuates into a qq 
pair, which scatters diffractively on the other proton. In such a reaction, only states with the same 
quantum numbers as the photon can be produced (the p, uj, 0, J and T mesons, their excitations 
and the Z boson). While vector meson photoproduction has been extensively studied in ep collisions 
at HERA, it has only recently been observed (for J/ij) and V'(2s) states) in hadron-hadron collisions 
at the Tevatron [145] (see section 7.2). It is not our aim here to present a complete survey of the 
physics which can be studied via exclusive photoproduction at the LHC. Instead, we focus on two 
examples, T production and Z production. 

At the LHC, photoproduced T — > M^A*~ "^^J be measurable, with one proton detected 
(there is no acceptance to measure both). The cross section is sensitive to QCD saturation effects, 
which arise as a result of non-linear gluon dynamics. The large mass of the T offers hope that 
this effect could be quantified in QCD perturbation theory. Central exclusive T production at 
the LHC has been studied in Refs. [147, 146, 216, 217, 218, 219]. The study in Ref. [219] presented 
results using a model for diffractive scattering which fits both the existing HERA data, not only for T 
production but also for the total inclusive deep inelastic scattering cross section. The model includes 
a rudimentary modelling of saturation effects and predicts a cross section of 5.1 pb for 14 TeV pp 
collisions at the LHC assuming neither proton is detected. This cross section is for T IJ-'^ 
with pf^ > 4 GeV/c. The cross section falls to 3.0 pb if one proton is detected at 420 m from the 
interaction point. These are large cross sections and can be measured with early LHC data. Tagging 
one proton should help control the pile-up background, and proton dissociation e.g. p — >■ pn^'K~ , 
thereby making the measurement possible also at high luminosities. The large rates mean that it will 
also be possible to measure the upsilon rapidity and this should help distinguish between different 
models of high-energy scattering. Measuring one proton makes possible a direct measurement of the 
photoproduction cross section, ct(7p -t- Tp) since a cut on the transverse momentum (e.g. around 
100 — 300 MeV) of the scattered proton will permit one to divide the event sample into two halves: 
one in which the detected proton predominantly radiated the photon and one in which it was the 
undetected proton that radiated the photon. The low pt sample corresponds to a measurement 
at higher values of the 7p centre-of-mass energy, with sensitivity up to 1 TeV in the 7p centre-of- 
mass energy. The prediction of the photoproduction cross section is, however, complicated by any 
dependence of gap survival on the scattered proton pT [220] . On the other hand, the measurement 

gives information on S"^ . It should also not be forgotten that central production of vector mesons 
can proceed also via odderon exchange, and it would be very interesting to observe that. 

Also possible at the LHC is the diffractive photoproduction of Z bosons [157, 147]. As 
discussed in section 7.2, the Standard Model cross section for Z production is negligibly small at 
the Tevatron: d(j/dy|y=o = 0.077 fb, compared with an upper limit from CDF of 0.25 pb at 95% 
C.L. [156]. At the LHC however, the cross section is expected to be about 6 fb for |y(^)| < 2, 
or 40 exclusive Z — )• e+e" or ji^ ix~ events per 100 fb~^ [147]. Clearly one will have to work at 
high luminosities, which means that pile-up backgrounds need bringing under control. This should 
be possible by selecting events with only the Z+Z" tracks on the vertex, a dilcpton invariant mass 
matching that of the Z, small pt{1'^1~) and tt — A0(Z+Z~), in conjunction with the requirement that 
the vertex be coincident with that obtained using two measured protons (at 420 m) . 
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9 Concluding remarks 



The subject of central exclusive production promises to deliver an entirely new dimension to the 
physics that will be explored at the LHC. With the installation of forward low-angle proton detectors, 
central masses can be measured out to a few hundred GeV/c^ and to a precision at or below 

1 GcV/c^, even in the presence of missing energy, using very few events. Moreover, tagging the 
protons effectively filters the central system with a strong bias towards J = 0, even-parity systems. 
Such a filter will provide a unique handle on the nature of any new physics that might be produced. 
It is true to say that this programme constitutes a huge advance over previous studies of CEP 
performed mainly at the ISR, SppS and Tevatron. For the first time (with the exception of the 
glueball searches), central masses can be produced in a range where interesting new phenomena are 
expected. 

Measurements from CDF at the Tevatron are very encouraging and support the validity 
of the theoretical calculations based in perturbative QCD: The theory is probably not too far off 
the mark. A second highlight of the recent past has been the demonstration that high luminosity 
backgrounds can be brought under control, which opens up the possibility to explore low (i.e. below 
1 fb) cross section processes at the LHC. 

There are many concrete examples of new physics that may reveal itself in CEP, and we 
have reviewed many of them here, but it should always be understood that CEP is principally a 
direct source of glue-glue (or photon-photon) interactions with a filter on the final state quantum 
numbers. It is thus ideally suited to a study of any new physics that couples to gluons (or photons). 
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